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Abstract. The billiard motion inside an ellipsoid Q C M"+^ is completely integrable. Its 
phase space is a symplectic manifold of dimension 2n, which is mostly foliated with Liouville 
tori of dimension n. The motion on each Liouville torus becomes just a parallel translation 
with some frequency that varies with the torus. Besides, any billiard trajectory inside Q is 
tangent to n caustics , • ■ • , Qa„, so the caustic parameters A = (Ai, . . . , A„) are integrals 
of the billiard map. The frequency map A w is a key tool to understand the structure of 
periodic billiard trajectories. In principle, it is well-defined only for nonsingular values of the 
caustic parameters. 

We present two conjectures, fully supported by numerical experiments. We obtain, from 
one of the conjectures, some lower bounds on the periods. These bounds only depend on the 
type of the n caustics. We describe the geometric meaning, domain, and range of w. The 
map u) can be continuously extended to singular values of the caustic parameters, although it 
becomes "exponentially sharp" at some of them. 

Finally, we study triaxial ellipsoids of M"^. We compute numerically the bifurcation curves 
in the parameter space on which the Liouville tori with a fixed frequency disappear We 
determine which elUpsoids have more periodic trajectories. We check that the previous lower 
bounds on the periods are optimal, by displaying periodic trajectories with periods four, five, 
and six whose caustics have the right types. We also give some new insights for ellipses of M?. 
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1. Introduction 

Birkhoff [7] introduced the problem of convex billiard tables more than 80 years ago as a way 
to describe the motion of a free particle inside a closed convex curve such that it reflects at 
the boundary according to the law "angle of incidence equals angle of reflection". He also 
realized that this billiard motion can be modeled by an area preserving map defined on an 
annulus. There exists a tight relation between the invariant curves of this billiard map and 
the caustics of the billiard trajectories. Caustics are curves with the property that a trajectory, 
once tangent to it, stays tangent after every reflection. Good starting points in the literature of 
the billiard problem are [|28ll39ll26ll . We also refer to [27] for some nice figures of caustics. 
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When the billiard curve is an ellipse, any billiard trajectory has a caustic. The caustics 
are the conies confocal to the original ellipse: confocal ellipses, confocal hyperbolas, and the 
foci. The foci are the singular elements of the family of confocal conies. In this case, the 
billiard map is integrable in the sense of Liouville, so the annulus is foliated by its invariant 
curves, the billiard map becomes just a rigid rotation on its regular invariant curves, and the 
rotation number varies analytically with the curve. 

The billiard dynamics inside an ellipse is known. We stress just three key results related 
with the search of periodic trajectories. First, Poncelet showed that if a billiard trajectory 
is periodic, then all the trajectories sharing its caustic conic are also periodic [|33ll . Second, 
Cayley gave some algebraic conditions to determine the caustic conies whose trajectories 
are periodic iH. Third, the rotation number can be expressed as a quotient of elliptic 
integrals [1311 |40l |43l . We note that the search of periodic trajectories inside an ellipse can 
be reduced to the search of rational rotation numbers. 

A rather natural generalization of this problem is to consider the motion of the particle 
inside an ellipsoid of ]R"+^. Then the phase space is no longer an annulus, but a symplectic 
manifold of dimension 2n. Many of the previous results have been extended to ellipsoids, 
although those extensions are far from being trivial. For instance, any billiard trajectory inside 
an ellipsoid has n caustics, which are quadrics confocal to the original ellipsoid. This situation 
is fairly exceptional, since quadrics are the only smooth hypersurfaces of ]R"+^, n > 2, 
that can have caustics flU. Then the billiard map is still completely integrable in the sense 
of Liouville, being the caustics a geometric manifestation of its integrability. In particular, 
the phase space is mostly foliated with Liouville tori of dimension n. The motion on each 
Liouville torus becomes just a parallel translation with some frequency that varies with the 
torus. Some extensions of the Poncelet theorem can be found in [[HI [121 [IS IMl. Several 
generalized Cay ley-like conditions were stated in lfT5l[T6l[T7l[T8ll . Finally, the frequency map 
can be expressed in terms of hyperelliptic integrals, see [fT3ll34l . The setup of these last two 
works is M?, but their formulae are effortless extended to ]R"+^. 

From Jacobi and Darboux it is known that hyperelliptic functions play a role in the 
description of the billiard motion inside ellipsoids and the geodesic flow on ellipsoids. 
Nevertheless, we skip the algebro-geometric approach (the interested reader is referred 
to [|32ll29l[30l l2ll3l) along this paper, in order to emphasize the dynamical point of view. Here, 
we just mention that the billiard dynamics inside an ellipsoid can be expressed in terms of 
some Riemann theta-functions associated to a hyperelliptic curve, and so, one can write down 
explicitly the parameterizations of the invariant tori that foliate the phase space; see pni22l| . 

Periodic orbits are the most distinctive special class of orbits. Therefore, the first task to 
carry out in any dynamical system should be their study, and one of the simplest questions 
about them is to look for minimal periods. In the framework of smooth convex billiards the 
minimal period is always two. Nevertheless, since all the two-periodic billiard trajectories 
inside ellipsoids are singular — in the sense that some of their caustics are singular elements 
of the family of confocal quadrics — , two questions arise. Which is the minimal period among 
nonsingular billiard trajectories? Which ellipsoids display such trajectories? 

In order to get a flavor of the kind of results obtained in this paper, let us consider the 
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three-dimensional problem. Let Q be the triaxial ellipsoid given by /a + + /c = 1, 
with < c < b < a. We assume that a = 1 without loss of generality. Any billiard trajectory 
inside Q has as caustics two elements of the family of confocal quadrics given by 

r x'^ xp' 

Qx = { {x,y,z) G : + + = 1 

[ a — A — A c — A 

We restrict our attention to nonsingular trajectories. That is, trajectories whose caustics are 
ellipsoids: < A < c, 1 -sheet hyperboloids: c < A < 6, or 2-sheet hyperboloids: 6 < A < a. 
The singular values A G {a, b, c} are discarded. It is known that there are only four types 
of couples of nonsingular caustics: EHl, HlHl, EH2, and H1H2. The notation is self- 
explanatory. It is also known that any nonsingular periodic billiard trajectory inside Q has 
three so-called winding numbers uiq, mi, m2 G N which describe how the trajectory folds in 
M^. For instance, mo is the period. The geometric meanings of mi and m2 depend on the type 
of the couple of caustics, see section [5| We stand out two key observations about winding 
numbers. First, some of them must be even. Namely, the ones that can be interpreted as the 
number of crossings with some coordinate plane. Second, we conjecture that they are ordered 
as follows: m2 < mi < ttiq. This unexpected behaviour is supported by extensive numerical 
experiments. In fact, we believe that it holds in any dimension. The combination of both 
observations crystallizes in the following lower bounds. 

Theorem 1. If the previous conjecture on the winding numbers holds, any periodic billiard 
trajectory inside a triaxial ellipsoid ofM.^ whose caustics are of type EHl, HI HI, EHl, and 
HlHl has period at least five, four, five, and six, respectively. 

All the billiard trajectories of periods two and three are singular. The two-periodic ones 
are contained in some coordinate axis, so they have two singular caustics. The three-periodic 
ones are contained in some coordinate plane, so they have one singular caustic. 

We shall prove in section [3] the generalization of these lower bounds to any dimension, 
see theorem |9j Samples of periodic trajectories with minimal periods are shown in figure 13 



Hence, these lower bounds are optimal. Next, we look for ellipsoids with minimal periodic 
trajectories. We recall that a = 1, so each ellipsoid Q is represented by a point in the triangle 
P = {(b,c) G : < c < 6 < 1}. Let Pi*, P^, P^, and P^ be the four regions of P 
that correspond to ellipsoids with minimal periodic trajectories of type EHl, HlHl, EH2, and 
H1H2, respectively. They are shown in figure[T| Their shapes are described below. 

Numerical Result 1. Letr = (3 - v^)/2 ^ 0.382, 6* = b^ = 1, andbl = bl = 1/2. Then 

P; = {{b,c)eP:b<b*, c<g*{b)}, 1<J<4, 
for some continuous functions g* : [0, b*] — M such that 

(i) gl is concave increasing in [0, 1], < gl{b) < bfor all b G (0, 1), and gl{l) = r; 

(ii) (?2 concave increasing in [0, 1], gl{b) < g2{b) < bfor all b G (0, 1), and g^^iX) = 1/2.' 

(Hi) gl is the identity in [0, r], concave decreasing in [r, 1/2], and gl{l/2) = 0; and 

(iv) gl is increasing in [0, 1/3], concave decreasing in [1/3, 1/2], 36/4 < gl{b) < bfor all 
b G (0,1/3), < glib) < g*{b) for all b e (1/3,1/2), ^^(1/3) = 1/4, and gl{l/2) = 0. 




b b 

Figure 1. The four regions of ellipsoids with minimal periodic txajectories. Left: The yellow 
region (type EHl, period 5) is contained in the green one (type HlHl, period 4). Right: The 
blue region (type H1H2, period 6) is contained in the magenta one (type EH2, period 5). 



The functions g* can be explicitly expressed by means of algebraic formulae. We shall 



prove that 5^2 (^) = 6/(1 + &) in proposition 18 We shall study the other three functions 



glib) 



in another paper BTl . because the techniques change drastically. A generalized Cay ley-like 
condition is the main tool. For instance, 

(1 - 6/2 - y/b{l -3b/A))b/{l - bf for < 6 < 1/3 
1-26)6/(1-6)2 for 1/3 < 6 < 1/2 ■ 

This function g*^ is not concave in the interval [0,1/3]. 

We shall describe in section [5] the regions corresponding to ellipsoids that have periodic 
trajectories with given winding numbers (or quasiperiodic trajectories with given frequencies) 
for the four caustic types. Those general regions have the same shape as these "minimal" 
regions. That is, they are below the graphs of some functions with properties similar to 
the ones stated previously. Therefore, we discover a general principle. The more spheric 
is an ellipsoid, the poorer are its four types of billiard dynamics. Here, spheric means 



(6, c) ~ (1, 1). We quantify this principle in propositions 15 and 17 



The key step for the numerical computation of these regions is to explicitly extend the 
frequency map for singular values of the caustic parameters. The extension is "exponentially 
sharp" at some points, which implies another general principle. The billiard trajectories with 
some almost singular caustic are ubiquitous. We shall enlighten it in subsection |5.7| by giving 
a quantitative sample. The minimal periodic trajectories shown in figure 13 reinforce it. 

Finally, we want to mention that there exist many remarkable results about periodic 
trajectories in other billiard and geodesic problems. For instance, several nice algebraic 
closed geodesies on a triaxial ellipsoid can be seen in [|23l . and a Cay ley-like condition for 
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billiards on quadrics was established in [1]. Some results stray from any integrable setup. 
For example, some general lower bounds on the number of periodic billiard trajectories inside 
strictly convex smooth hypersurfaces can be found in [|4l|20l|2T]|. The planar case was already 
solved by Birkhoff [7J. Of course, these lower bounds are useless for integrable systems 
where the periodic trajectories are organized in continuous families. 

We complete this introduction with a note on the organization of the paper. In section [2] 
we review briefly some well-known results about billiards inside ellipsoids in order to fix 
notations that will be used along the rest of the paper. Next, the frequency map is introduced 
and interpreted in section [s} This section, concerning ellipsoids of ]R"+^, also contains two 
conjectures and the lower bounds on the periods. Billiards inside ellipses of and inside 
triaxial ellipsoids of are thoroughly studied in sections |4] and [sj respectively. Billiards 
inside nondegenerate ellipsoids of ]R"+^ are revisited in section [6| Some technical lemmas 
have been relegated to the appendices. 



2. Preliminaries 



In this section details are scarce and technicalities are avoided. Experts can simply browse 
this section. We will list several basic references for the more novice readers. 



2.1. Confocal quadrics and elliptic billiards 

The following results go back to Jacobi, Chasles, Poncelet, and Darboux. 

The starting point of our discussion is the n-dimensional nondegenerate ellipsoid 

{n+l 2 1 
a; = (xi, . . . , Xn+i) e : E ^ = 1 1 ' 

where ai, . . . , a„+i are some fixed real constants such that < ai < ■ ■ ■ < a^+i. The 
degenerate cases in which the ellipsoid has some symmetry of revolution are not considered 
here. This ellipsoid is an element of the family of confocal quadrics given by 

f ^ t2 1 

The meaning of is unclear in the singular cases jj, G {ai, . . . , a„+i}. In fact, there are two 
natural choices for the singular confocal quadric when fx = aj. The first choice is to define 
it as the n-dimensional coordinate hyperplane 

Hj = {x = (xi, . . . , Xn+i) G : x, = 0} , 

but it also makes sense to define it as the {n — 1) -dimensional /oca/ quadric 

Fj = \x = {xi,..., Xn+i) G ■.Xj = and ^ — = 1 

which is contained in the hyperplane Hj . Both choices fit in the framework of elliptic billiards, 
but we shall use the notation Qa^ = Hj along this paper. 
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Theorem 2 ( [i32ll28l l2ll39ll). Once fixed a nondegenerate ellipsoid Q, a generic line i C 

is tangent to exactly n distinct nonsingular confocal quadrics Qx^, • • • , Qx„ such that Ai < 

A2 < ■ ■ • < A„, Ai G (-cx),ai) U (01,02), and \i G (aj_i,ai) U {at, ai+i), for i = 2, . . . ,n. 

Set ao = 0. If a generic line i has a transverse intersection with the ellipsoid Q, then 
Ai > 0, so Ai G (flO) c^i) U (ai, 02). The value Ai = is attained just when i is tangent to Q. 
A line is generic in the sense of the theorem if and only if it is neither tangent to a singular 
confocal quadric[|]nor contained in a nonsingular confocal quadric. 

If two lines obey the reflection law at a point x G Q, then both lines are tangent to 
the same confocal quadrics f39^. This shows a tight relation between elliptic billiards and 
confocal quadrics: all lines of a billiard trajectory inside the ellipsoid Q are tangent to exactly 
n confocal quadrics Qx^, ■ ■ ■ , Qx,,, which are called caustics of the trajectory. We will say 
that A = (Ai, . . . , A„) G M" are the caustic parameters of the trajectory. 

Definition 1. A billiard trajectory inside a nondegenerate ellipsoid of the Euclidean space 
]R"+^ is nonsingular when it has n distinct nonsingular caustics; that is, when its caustic 
parameter belongs to the nonsingular caustic space 

A=|(Ai,...,A„)GM- °<//<^^<;V<^" . ^^.^ 1.(2) 
y Xi G (aj_i, ai) U [tti, ai+i)for l<i<n J 

We will only deal with nonsingular billiard trajectories along this paper. We denote the 
2" open connected components of the nonsingular caustic space as follows: 

A. = ((Ai,...,A„) G W : ° < ^ ^ . ^ 

y Ai G (fli+CT.-i, cti+aj for 1 < z < n 

for 0" = ((Ti, . . . , cr„) G {0, 1}"^. For instance, the first caustic Qai is an ellipsoid if and only if 
Ai G (flO) cti); that is, if and only if A G A^. with ai = 0. We will draw the space A for ellipses 
and triaxial ellipsoids of in sections |4] and [5} respectively. 

Theorem 3. If a nonsingular billiard trajectory closes after mo bounces, all trajectories 
sharing the same caustics close after itlq bounces. 

Poncelet proved this theorem for conies [33] . Darboux generalized it to triaxial ellipsoids 
of M'^. Later on, this result was generalized to any dimension in [iTTl[T2l[T3ll35l . 

Theorem 4. The nonsingular billiard trajectories sharing the caustics Qx^ , ■ ■ ■ , Qa„ close 
after m bounces — up to the action of the group of symmetries G = (Z/2Z)'^+^ ofQ — , if and 
only ifm > n + 1 and 

^ hm+l ■ ■ ■ hn+2 ^ 

rank '■ '■ < m — n, (3) 

\ ^2m-l ■ ■ ■ hn+m ) 

where ^/(ai - s) ■ ■ ■ (a„+i - s)(Ai - s) ■ ■ • (A„ - s) = /iq + + ^2S^ H • 

% By abuse of notation, it is said that a line is tangent to the singular confocal quadric Qa^ when it is contained 
in the coordinate hyperplane Hj or when it passes through the focal quadric Fj. 
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The group G is formed by the 2"+^ reflections — involutive linear transformations — with 
regard to coordinate subspaces. The phrase "a billiard trajectory closes after m bounces up 
to the action of G" means that if {qk)k& is the sequence of impact points of the trajectory, 
then there exists a reflection g E G such that qk+m = g{qk) for all k E Z. Hence, billiard 
trajectories that close after m bounces up to the action of the group G, close after mg = m or 
mo = 2m bounces, because qk+2m = g{qk+m) = g^{qk) = qk- 

Cay ley proved this theorem for conies ||9l. Later on, this result was generahzed to any 
dimension by Dragovic and Radnovic in lfT5lfT6l . 

2.2. Complete integrability of elliptic billiards 

We recall some results obtained by Liouville, Arnold, Moser, and Knorrer. 

A symplectic map f : M ^ M defined on a 2n-dimensional symplectic manifold is 
completely integrable if there exist some smooth /-invariant functions Ji, . . . , J„ : M — )■ M 
(the integrals) that are in involution — that is, whose pair-wise Poisson brackets vanish — and 
that are functionally independent almost everywhere on the phase space M. In this context, 
the map / = (Ji, . . . , /„) : M — )• is called the momentum map. A point m G M is a 
regular point of the momentum map when the n-form d/i A ... A d/„ does not vanish at m. 
A vector A G is a regular value of the momentum map when every point in the level set 
/^^(A) is regular, in which case the level set is a Lagrangian submanifold of M and we say 
that /~^(A) is a regular level set. 

The following result is a discrete version of the Liouville- Arnold theorem. 

Theorem 5 ([[421). ^ny compact connected component of a regular level set /^^(A) is 
diffeomorphic to T", where T = M/Z. In appropriate coordinates the restrictions of the 
map to this torus becomes a parallel translation 6 ^ 6 + uj. The map X u is smooth at the 
regular values of the momentum map. 

Thus the phase space M is almost foliated by Lagrangian invariant tori and the map on 
each torus is simply a parallel translation. These tori are called Liouville tori, the shift to is the 
frequency of the torus, and the map A i— )■ w is the frequency map. The dynamics on a Liouville 
torus with frequency to is mg-periodic if and only if niQUJ E Z". Liouville tori become just 
invariant curves when n = 1, in which case the shift is usually called the rotation number of 
the invariant curve, and denoted by p, instead of cu. 

Now, let Q be a (strictly) convex smooth hypersurface of ]R"+^ diffeomorphic to the 
sphere not necessarily an ellipsoid. The billiard motion inside Q can be modelled by 
means of a symplectic diffeomorphism defined on the 2ri-dimensional phase space 

M = {{q, p) E Q X S"' : p is directed outward Q atq} . (4) 

We define the billiard map f : M M, f{q,p) = iq',p'), as follows. The new velocity 
p' is the reflection of the old velocity p with respect to the tangent plane TgQ. That is, if we 
decompose the old velocity as the sum of its tangent and normal components at the surface: 
p = Pt + Pn with pt E TqQ and pn E NgQ, then p' = Pt — Pn = P — '^Pn- The new impact 
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point q' is the intersection of the ray {g + np' : /i > 0} with the surface Q. This intersection 
is unique and transverse by convexity. 

Elliptic billiards fit in the frame of the Liouville- Arnold theorem. 

Theorem 6 ( [i32lf30ll4Tl l2l). The billiard map associated to the nondegenerate ellipsoid ([7]) 
is completely integrable and the caustic parameters Ai, . . . , A„ are the integrals. The set of 
regular values of the corresponding momentum map is given by (|2]). 



3. The frequency map 

3.1. Definition and interpretation 

The rotation number for the billiard inside an ellipse is a quotient of elliptic integrals; 
see ll3Tl[TT1l . Explicit formulae for the frequency map of the billiard inside a triaxial ellipsoid 
of can be found in [fT3l §III.C]. An equivalent formula is given in [|34l §5]. Both formulae 
contain hyperelliptic integrals and they can be effortless generalized to any dimension. Since 
we want to avoid as many technicalities as possible, we will not talk about Riemann surfaces, 
basis of holomorphic differential forms, basis of homology groups, period matrices, or other 
objects that arise in the theory of algebraic curves. 

The following notations are required to define the frequency map. Once fixed the 
parameters ai, . . . , a„+i of the ellipsoid, and the caustic parameters Ai, . . . , A„, we set 

2n+l 

T{s) = ]^ (cj - s), {ci, . . . , C2„+i} = {ai, . . . , a„+i} U {Ai, . . . , A„}. 

i=l 

If A G A, then ci, . . . , C2n+i are pair-wise distinct and positive, so we can assume that 

Co := < Ci < ■ ■ ■ < C2„+i. (5) 
Hence, T(s) is positive in the n + 1 open intervals (c2j, C2j+i), and the improper integrals 

Kij= I ^ , i = 0,...,n-l, j = 0,...,n (6) 

are absolutely convergent, real, and positive. We also consider the n + 1 column vectors 

ir, = (iroi,...,i^n-i,,)*eM". 

It is known that vectors i^i, . . . , Kn are linearly independent; see [|25l §111.3]. 

Definition 2. The frequency map : A — )■ M" o/ the billiard inside the nondegenerate 
ellipsoid Q is the unique solution of the system ofn linear equations 

n 

K^ + 2^{-iyujjKj = {). (7) 

i=i 

Remark 1. Sometimes it is useful to think that the frequency u depends on the parameter 
c = (ci, . . . , C2n+i) G ]R+"^^, and not only on the caustic parameter A = (Ai, . . . , A„) G A. In 
such situations, we will write u = zu{c). This map c i— i- zu(c) is homogeneous of degree zero 
and analytic in the domain defined by inequalities Q. Homogeneity is deduced by performing 
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a change of scale in the integrals ([6]). Hence, we can assume without loss of generality that 
C2n+i = CLn+i = 1- Aualyticity follows from the fact that the integrands in (|6]) are analytic 
with respect to the variable of integration in all the intervals of integration and with respect to 
c, as long as condition ([5]) takes place. 

This definition coincides with the formulae contained in [[T3l l34l for n = 2. It 
is motivated by the characterization of periodic billiard trajectories contained in the next 
theorem. The factor 2 has been added to simplify the interpretation of the components of 
the frequency map, which are all positive, due to the factors (—1)^ . 

Theorem 7 ( ifTTl [HI). The nonsingular billiard trajectories inside the nondegenerate 
ellipsoid Q are periodic with exactly rrij points at Qc^^ and nij points at Qcj^+i if and only if 

moi^o+E;=i(-i)''^.^j = o. 

The numbers uiq, mi, ... , m„ that appear in theorem |7] are called winding numbers. The 
nonsingular billiard trajectories with caustic parameter A are periodic with winding numbers 

mo, mi, . . . , m„ if and only if 

in ■ 

^.■W = 7r^eQ+, j = i,...,n. (8) 

2mo 

We note that mo is the number of bounces in the ellipsoid Q = Qcq, so it is the period. 
Remark 2. The sequence of winding numbers of a nonsingular periodic billiard trajectory 
contains information about how the trajectory folds in the space IR"^"^ The following 
properties can be deduced from results contained in [ , 17. §4.1]. Here, "number of _" means 
"number of times that any periodic billiard trajectory with those caustic parameters do _ along 
one period". The intervals (c2j, C2j+i) with j 7^ can adopt exactly four different forms, each 
one giving rise to its own geometric picture. 

(i) If (c2j, C2j+i) = (ttj, then nij is the number of crossings with Hj, so it is even and 
mj/2 is the number of oscillations around the hyperplane Hj; 

(ii) If (c2j, C2j+i) = (Aj, ctj+i), then nij is the number of crossings with i/j+i, so it is even 
and mj/2 is the number of oscillations around the hyperplane Hj+i; 

(iii) If (c2j, C2j+i) = (oj, ttj+i), then nij is the number of (alternate) crossings with Hj and 
ifj+i, so it is even and nij/2 is the number of rotations that the trajectory performs when 
projected onto the {xj, x^+i) -coordinate plane tt^; and 

(iv) If (c2j, C2j+i) = {Xj, Aj+i), then m^ is the number of (alternate) tangential touches with 
Qxj and Qa^+i , so it can be even or odd, and it is the number of oscillations between both 
caustics. 

These four properties suggest us the following definitions, which establish the geometric 
meaning of the components of the frequency map. They change with the open connected 
components of the nonsingular caustic space. 

Definition 3. Let u = {ui, . . . , w^) : A — M" Z^e the frequency map. 

(i) If{c2j, C2j+i) = {dj, Aj+i), then Uj = mj/2mQ is the i/j -oscillation number; 

(ii) If{c2j, C2j+i) = {Xj, Oj+i), then Uj = mj/2mQ is the /fj+i-oscillation number; 
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(Hi) If{c2j, C2j+i) = (ctj, dj+i), then uj = mj/2mQ is the ttj -rotation number; and 
(iv) If{c2j, C2j+i) = {Xj, 2ujj = nij/mo is the {Qxj, Qa^+i) -oscillation number. 

Remark 3. It is important to notice that (only) when a mo-periodic billiard trajectory has two 
caustics of the same type — that is, when some interval (c2j, C2j+i) falls into the fourth case — , 
it is possible that mocu ^ Z", although then 2moUJ E Z". This is due to the factor 2 that we 
have added in the definition of the frequency map. 

Finally, if Qxj is not an ellipsoid — that is, if Ai > ai — , then ci = Oi, and tuq is the 
number of crossings with Hi, so it is even. Therefore, the following corollary holds. 

Corollary 8. Among all the nonsingular billiard trajectories inside a nondegenerate ellipsoid, 
only those with an ellipsoid as caustic can have odd period. 

3.2. Two conjectures 

Conjecture 1. The frequency map is a local diffeomorphism; i.e., it is nondegenerate: 



This conjecture has several relevant consequences throughout the paper. Popov and 
Topalov ||34| have shown that the frequency map is almost everywhere nondegenerate when 
Q is a triaxial ellipsoid of M?, although they only consider the components A^r such that 
cTi = 0. The nondegeneracy of the frequency map is important because it is an essential 
hypothesis — although we acknowledge that it can be replaced by some weaker Riissmann- 
hke nondegeneracy conditions [l38l §2] — in most KAM-like theorems, which are the standard 
tool to prove the persistence of Liouville tori under small smooth perturbations of completely 
integrable maps. Therefore, if conjecture [T] holds, we can conclude that most of the Liouville 
tori of the billiard phase space persist under small smooth perturbations of the ellipsoid. We 
shall present evidences for this conjecture in sections]?] and [5| 

Conjecture 2. Winding numbers of nonsingular periodic billiard trajectories are ordered in 
a strict decreasing way. More concretely. 



Inequality m„ > 2 is immediate, because C2„+i = ctn+i: so m„ is even. Inequalities 
rrij < nio for j > 1 are also immediate, because the number of crossings with any fixed 
hyperplane or the number of tangential touches with any fixed caustic can not exceed the 
number of segments of the periodic billiard trajectory. The strict inequalities mj < tuq could 
be also established (using the symmetries of the ellipsoid), but we skip the details, since such 
a small improvement is not worth the effort. 

If both conjectures hold, the components of the frequency map are also ordered in a strict 
decreasing way and lie in the range (0, 1/2); that is. 




VA e A. 



2 < rrin < ■ ■ ■ < nil < = period. 




To prove this, we note that conjecture[T]implies that ui{A) is open in M" and lo ^(Q") is dense 
in A, whereas conjecture ]2] and relation (]8]) imply that the strict inequalities < w„ < ■ ■ ■ < 
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ui < 1/2 hold for rational frequencies. Therefore, < c<;„(A) < ■■■ < wi(A) < 1/2 for 
any A G A, but these inequalities must be strict because a; (A) is open. We have numerically 
checked that inequalities (|9]) hold for thousands of random choices of ai, . . . , Ai, . . . , An 
in "dimensions" n <h. The details of the experiments for n = 2 are presented in section [5] 

3.3. Lower bounds on the periods 

We know from theorem |4] that the period niQ of any nonsingular periodic billiard trajectory 
inside a nondegenerate ellipsoid Q C ]R"+^ verifies that mo > n + 1. This result can be 
improved in several ways using the ordering of the winding numbers stated in conjecture |2] 
For instance, the global lower bound mo > n + 2 follows directly. We present below more 
refined semi-global lower bounds, holding each one on a different open connected component 
of the nonsingular caustic space. They are obtained by realizing that some winding numbers 
must be even in agreement with the first items in remark |2j The lower bound associated to 
some connected component reaches the value 2n + 2, which doubles the original lower bound 
given in theorem |4j 

Theorem 9. Given any a = (ai, . . . , cr„) G {0, 1}", let E^r be the subset o/{0, 1, . . . , n} such 
that: 1) e E„ ^ ai = 1; 2) j e ^ (0-^,0-^+1) 7^ (1, 0); and 3)n e E„. 

(i) If mo, . . . ,mn are the winding numbers of a periodic trajectory with caustic parameter 
in Ag-, then rrij is even for all j G E„. 

(ii) If conjecture^holds, any periodic billiard trajectory inside a nondegenerate ellipsoid of 
K'^+i whose caustic parameter is in K„ has period at least 

, , .1 3 2 < nin < ■ • ■ < mo sequence of integers 
x((t) := mm < mo : , , . ^ 

I such that nij is even for any j G E^ 

(Hi) Le? 1 = (1, . . . , 1) G {0, 1}", ^ = (..., 0, 1, 0, 1, 0) G {0, 1}", and a G {0, 1}". Then 
n + 2 = x(<;) < x(a) < x(l) = 2n + 2, ^a^l,q. 

Proof, (i) We recall that nij must be even in the three first cases listed in remark [2j This is 
the key property. For instance, m„ is always even because C2n+i = o-n+i- If = 1^ then 
Ai G (01,02) and ci = ai, so mo is even. If nij is odd, then {c2j,C2j+i) = {Xj,Xj+i) and 
Xj, Xj+i G {aj, flj+i), so (cTj, cTj+i) = (1, 0). Hence, we have seen that (aj, cTj+i) 7^ (1, 0) =^ 
rrij is even. 

(ii) This follows directly from the previous item and the definition of x(cr). 

(iii) First, we note that E^ = {0, . . . , n} and E^ = {. . . , n — 4, n — 2, n}. Therefore, 
x(l) = min {mo : 3 2 < m„ < • ■ ■ < mo sequence of even numbers} = 2n + 2, and 

, , . I 3 2 < m„ < ■ ■ ■ < mo sequence of 1 
x(<;j = mm < mo : > = n + 2. 

I mtegers s.t. m„, m„_25 • • • are even j 

The minimum value of mo among all integer sequences such that 2 < m„ < ■ • • < mo is 
attained at the sequence mj = n + 2 — j,0<j<n. Thus, x(cr) > n + 2 for all a G {0, 1}", 
and x(cr) = n + 2^E^ = E^^a = <^. 
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On the other hand, E^r C E„i =^ >c{a) < K{a'). Hence, x((t) < x(l) = 2n + 2 for all 
a e {0, 1}". Finally, = 2n + 2 ^ E„ = {0, . . . ,n} ^ aj Vj ^ a = 1. □ 

Remark 4. Let = (0, . . . , 0) and ^ = (...,1,0,1,0,1). Then = {l,...,n}and 
E^ = {. . . , n — 5, n — 3, n — 1, n}, so x(0) = 2n + 1 and x(f) = n + 3. 

Remark 5. All the lower bounds x((j) can be exphcitly computed when = 2. In that case, 

4 = x(0 < x(0) = 5 = x(^) < x(l) = 6 

where <; = (1, 0), = (0, 0), f = (0, 1), and 1 = (1, 1). Therefore, theorem [T] about triaxial 
ellipsoids of is just a particular case of theorem |9j It suffices to realize that a = 0, a = <j, 
(J = and (7 = 1 correspond to the cases EHl, HlHl, EH2, and H1H2, respectively. 

Remark 6. The function x : {0, 1}" {n + 2, . . . , 2n + 2} is surjective and has average 

x:=2"" x((t) =3n/2 + 2. 

o-e{o,i}" 

We skip the details, the proof is by induction over n. Thus, these semi-global lower bounds 
improve the global lower bound n + 2 by, in average, almost 50%. 

Now, a natural question arises. Are these semi-global lower bounds optimal? Optimal 
does not mean that there exists a x((t) -periodic billiard trajectory whose caustic parameter is 
in Ao- inside all nondegenerate ellipsoids, but just inside some of them. And we put another 
question. Which are the ellipsoids with such "minimal" periodic billiard trajectories? Both 
questions become almost trivial for ellipses; see subsection 4.6 The case of triaxial ellipsoids 
of was numerically answered in the introduction. The general case remains open, but we 
conjecture that all these semi-global lower bounds are optimal. 



4. Billiard inside an ellipse 

In this section we describe the main properties of the frequency map when n = 1, in which 
case it is called rotation number and denoted by p. Many of these properties are old, but 
the observation that the the rotation number is exponentially sharp at the singular caustic 
parameter seems to be new. The known results can be found in the monographes [[28ll39l and 
the papers [ISllIIIlllQllll. 



4.1. Confocal caustics 

To simplify the exposition, we write the ellipse as 

Q = hx,y) eR^ + j = l\ , a>b>0, 

where we could assume, without loss of generality, that a = 1; see remark [T] Then any 
nonsingular billiard trajectory inside Q is tangent to one confocal caustic of the form 
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where the caustic parameter A belongs to the nonsingular caustic spacej|] 

A = EUH, E = {0,b), H = {b,a). (10) 

We have chosen those names for the connected components of A because then Qx is an ellipse 
for A G -E, and a hyperbola for X e H. 

4.2. Phase portrait 

We describe now the billiard dynamics inside an ellipse. This description goes back to 
Birkhoff [7, §VIIL12], so it is rather old and we just list the results. Concretely, we want 
to know how the phase space is foliated by Liouville tori (invariant curves on which the 
motion becomes a rigid rotation) and separatrices (invariant curves on which the motion tends 
to some hyperbolic periodic trajectories). 

Let us put some global coordinates {if, r) over the billiard phase space M defined 
in (|4]), just for visualization purposes. First, following Birkhoff, we parameterize the impact 
points on the ellipse by means of an angular coordinate E T. We take, for instance, 
q = 7(^5) = (a^/^ cos V9, 6^/^ sin </)). Second, given an outward unitary velocity p G §, we 
set r = (7'((y9),p), and so |r| < |7'((y3)| = {asin^ ip + bcos^ cpY^'^. Then the correspondence 
(g, p) I— (if, r) allows us to identify the phase space M with the annulus 

A = {{ip,r) G T X M : < asinV + ^cosV} • (H) 

In these coordinates, the caustic parameter becomes \{lp, r) = (a — b) sin^ y9 + 6 — r^. The 
partition of the annulus into invariant level curves of A is shown in figure |2j 

Each regular level set contains two Liouville curves and represents the family of tangent 
fines to a fixed nonsingular caustic Qx- If Qx is an ellipse, each Liouville curve has a one- 
to-one projection onto the coordinate and corresponds to rotations around Qx in opposite 
directions, so they are invariant under /. If Qx is a hyperbola, then each Liouville curve 
corresponds to the impacts on one of the two pieces of the ellipse between the branches of 
Qx, so they are exchanged under / and invariant under 

The singular level set r) G A : X(lp, r) = b} gives rise to the 00-shaped curve 

X-\b) = {((^,r) G A : r = ±(a - b)^/^sm^} , 

which corresponds to the family of lines through the foci. This singular level set has rotation 
number 1/2; see ||26l page 428]. The cross points on this singular level represent the two- 
periodic trajectory along the major axis of the ellipse, and the eigenvalues of the differential 
of the billiard map at them are positive but different from one: e^'''' with cosh^ h/2 = a/b 
and h > 0. On the contrary, the two-periodic trajectory along the minor axis correspond to 
the centers of the regions inside the 00-shaped curve, and the eigenvalues in that case are 
conjugate complex of modulus one: e^^'^^^ with cos^ n9 = b/a and < 9 < 1/2. Therefore, 
the major axis is a hyperbolic (unstable) two-periodic trajectory and the minor axis is an 

§ When A ^ 6~ (resp., A 6+) the caustic Qx flattens into the region of the x-axis enclosed by (resp., outside) 
the foci of the elHpse Q. When A — > a^, the caustic flattens into the whole y-axis. 



The frequency map for billiards inside ellipsoids 14 




Figure 2. Phase portrait of the bilHard map in ((ys, r) coordinates for a = 1 and 6 = 4/9. 
The dashed black lines enclose the phase space ( [TT] ). The black points are the hyperbolic two- 
periodic points corresponding to the oscillation along the major axis of the ellipse. The black 
curves are the separatrices of these hyperbolic points. The magenta points denote the elliptic 
two-periodic points corresponding to the oscillation along the minor axis of the ellipse. The 
magenta curves are the invariant curves whose rotation number coincides with the frequency 
of these elliptic points. The invariant curves with rotation numbers 1/6, 1/4 and 1/3 are 
depicted in blue, green and red, respectively. The red points label a three-periodic trajectory 
whose caustic is an ellipse. The green points label a four-periodic trajectory whose caustic is 
a hyperbola. 



elliptic (stable) one. These are the only two-periodic motions. The basic results about the 
stability of two-periodic billiard trajectories can be found in [|28li39l . 

4.3. Extension and range of the rotation number 

Let p(A) be the rotation number of the billiard trajectories inside the ellipse Q sharing the 
nonsingular caustic Q\. From definition [2] we get that the function p : E H ^ given 
by the quotients of elliptic integrals 

, ^ ^{X-s){b-s){a-s) h ^t(t~l){t~x) 

p{\)= p{\-h,a) = d. = ^ .1 dt ' (12) 

Jmax(fe,A) ^{A-s)(6-s){a-s) JO ^t{t-l)(t-x) 

where the parameters 1 < x < are given by x = {a — m) / {a — Wi) and p = a/{a — m), 
with m = min(6, A) and m = max(6, A). The second equality follows from the change of 
variables t = {a — s)/{a — rn). The second quotient already appears in [|13l . Other equivalent 
quotients of elliptic integrals were given in [13111431 . We have drawn the rotation function p(A) 
in figure |3j compare with [43, figure 2]. 
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Figure 3. The rotation function p{\) of the ellipse for a — \ and b = 4/9. Colours are taken 
from figure|2] The parameters Aj. fast approach b as tends to 1/2. 



Proposition 10. The function p : E U H ^ M. given in ([72]) has the following properties. 

( i) It is analytic in A = E U H. 

(ii) It can be continuously extended to the closed interval A = A U dA = [0, a] with 

p(0) = 0, p{b) = 1/2, p{a) = Q, 

where the limit value Q < g < 1/2 is defined by sin^ txq = b/a. 
( Hi) Let k'^ and be the positive constants given by 

ds 



cosh K = a/b. 



lb \/ s[s — b){a — s] 
The asymptotic behavior of p{X) at the singular parameters A G OA = {0, b, a} is: 

(a) p(A) = + 0(A3/2), A ^ 0+; 

(b) p(A) = 1/2 + K^/ log |6 - A| + O (1/ log^ \b - \\), as X ^ b; and 

(c) p(A) = g + 0(a — A), as \ ^ a^. 

(iv) Given any p° G {g, 1/2), let A° be the biggest parameter in E such that p(A° ) = p", and 
let Al be the smallest parameter in H such that p(Al) = p". Both parameters become 
exponentially close to the singular caustic parameter b when po tends to 1/2. In fact, 



Xl = b± 16(a - 6)e-''"/(i/2-P°) + O 



-2KS/(l/2-pO) 



), P°^(l/2)'. 



Proof, (i) follows from remark[T] The rest of the proof is postponed to Appendix A. 2 □ 



Remark 1. If conjecture [T] holds, then p'(A) is positive in E and negative in H, so p(A) maps 
diffeomorphically E onto (0, 1/2) and H onto {g, 1/2). In particular, the parameters A° and 
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A']_ are unique. The conjecture remains unproven, but we shall see in proposition 1 1 that p(A) 
is increasing in E, which suffices to check the unicity of A° . 

Remark 8. The limit rotation number q is related to the conjugate complex eigenvalues e'^'^'^^'' 
of the elliptic two-periodic orbit. Concretely, 6 + q = 1/2. Besides, g tends to zero 
when the ellipse flattens and tends to one half when the ellipse becomes circular. That is, 

\imb/a^o+ = 0, and limb/a^i- Q = 1/2. 

Definition 4. The continuous extension p : [0, a] — t- M z5 called the (extended) rotation 
function of the ellipse Q. 

4.4. Geometric meaning of the rotation number 

Let us assume that the billiard trajectories sharing some nonsingular caustic Qx are mo- 
periodic, so they describe polygons with niQ sides inscribed in the ellipse Q. Then, according 
to theorem|7} equation ([S]), and corollary [8} it turns out that p(A) = mi/2mo for some integers 
2 < mi < mo such that mi is always even whereas mo can be odd only when Q\ is an ellipse. 
Besides, from the geometric interpretation of the frequency map presented in section |3} we 
know that: 1) If Qx is an ellipse, the polygons are enclosed between Q and Q\, and make 
mi/2 turns around the origin; and 2) If Q\ is a hyperbola, they are contained in the region 
delimited by Q and the branches of Q\, and cross mi times the minor axis of the ellipse. 
These interpretations can be extended to nonperiodic trajectories. Concretely, 

lim Hk/k if A G -E, 



P(A) 



fc- 



\ Um Ik/k if X e H, 

^ fc-s>+oo 



where (respectively, /fc) is the number of turns around the origin (respectively, crossings of 
the minor axis) of the first k segments of a given billiard trajectory with caustic Qa- 

Proposition 11. The rotation function p(A) is increasing in E. 

Proof. Let 7 : T — ?■ Q be a fixed parameterization of the ellipse Q. Then the billiard dynamics 
inside Q associated to any caustic Qx, \ E E, induces a circle diffeomorphism : T — > T 
of rotation number p(A). Let < Ai < A2 < 6. The billiard trajectories sharing the small 
caustic Qx2 rotate faster than the ones sharing the big caustic Qx^, so Fx-^ < Fx^ for any 
two compatible lifts Fx^ of the circle diffeomorphisms fxy Then p(Ai) < p(A2); see [|26l 
Proposition 11.1.8]. Thus, p(A) is nondecreasing and, by analyticity, increasing. □ 
We have not proved that p(A) is decreasing in H because it is not easy to construct an 
ordered family of circle diffeomorphisms for caustic hyperbolas. 

4.5. Bifurcations in parameter space 

We want to determine all the ellipses Q = {x^ / a + y^/6 = 1}, < 6 < a, that have billiard 
trajectories with a prescribed rotation number p° G (0, 1/2) and with a prescribed type of 



The frequency map for billiards inside ellipsoids 



17 




Figure 4. Examples of symmetric nonsingular billiard trajectories with minimal periods for 
a = 1 and 6 = 4/9. Left: Period three and the caustic is an ellipse. Right: Period four and the 
caustic is a hyperbola. The continuous lines are reserved for the trajectories that correspond to 
the periodic orbits depicted in figure |2] 



caustics (ellipses or hyperbolas). We recall that the rotation function p(A) diffeomorphically 
maps E onto (0, 1/2), and H onto {q, 1/2). Therefore, p° G p{E) for all ellipses Q, whereas 

p° G p{H) ^ g < p° sin^ ng < sin^ np^ -^b < a sin^ 7rp°. (13) 

This shows that flat ellipses have more periodic trajectories than rounded ones. There exist 



similar results for triaxial ellipsoids of M . See, for instance, propositions 15 and 
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4.6. Examples of periodic trajectories with minimal periods 

The billiard map associated to an ellipse has no fixed points, its only two-periodic points 
correspond to the oscillations along the major or minor axis, and only the trajectories with an 
ellipse as caustic can have odd period. Therefore, the periodic trajectories with an ellipse as 
caustic have period at least three, whereas the ones with a hyperbola as caustic have period at 
least four. These lower bounds are optimal; see figure |4j To be more precise, we set 

X* = 3a& ^ ^ 

^ a + b + 2 -ba + b^' " a - b' 

We note that e E for aU < 6 < a, and G if for aU < 6 < a/2. The trajectories 
with caustic Qx*^ are three -periodic, the ones with caustic Qx*^ are four-periodic. The proof is 
an elementary exercise in Euclidean geometry. We leave it to the reader. Finally, we deduce 
from the geometric interpretation of the rotation number given before that p{X^) = 1/3 and 



p(A^) = 1/4. This second identity explains the restriction b < a/2; see ( 13). 
5. Billiard inside a triaxial ellipsoid of 



The previous section sets the basis of this one. Roughly speaking, we want to follow the same 
steps — extension of the frequency map and description of its range — in order to find the 
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Oblate 
ellipsoids: 

h = l 



Segments: 



Circles: 



(6,c) = (0,0) 



Flat ellipsoids: c — Q 



(5,c) = (l,0) 



Figure 5. The triangular parameter space P. 



same results — bifurcations in the parameter space and minimal periodic trajectories. But the 
study of ellipsoids is harder, which has two unavoidable consequences. First, statements and 
proofs of the analytical results are more cumbersome. Second, some results remain unproven, 
so we shall present numerical experiments and semi-analytical arguments as support. 

5.1. Confocal caustics 

The caustics of a billiard inside a triaxial ellipsoid are described in several places. The 
representation of the caustic space shown in figure |6] can also be found in [|30ll44l[T9l . 
We write the triaxial ellipsoid as 



We could assume, again without loss of generality, that a = 1. Then the parameter space of 
triaxial ellipsoids in M? can be represented as the triangle 



whose edges represent ellipsoids with a symmetry of revolution (oblate and prolate ones) or 




a > 6 > c> 0. 



P= {(6,c) gM^:0<c<6<1} 



(15) 



flat ellipsoids, as illustrated in figure [5] We shall write the statements of the main results for 
arbitrary values of a, but we shall take a = 1 in the pictures. 



From theorem |2| we know that any nonsingular bilhard trajectory inside the ellipsoid Q 
is tangent to two distinct nonsingular caustics of the confocal family 



r iP' 1 
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Figure 6. The nonsingular caustic space K = {Ex Hi) [J{Ex H2) U {Hi Fi) U {Hi x H2) 
with its border = G U i? U 5 U A°. 



The caustic Q\ is an ellipsoid for A e -E, a one-sheet hj^erboloid when A e Hi, and a 
two-sheet hyperboloid if A e -ff2> where 

£; = (0,c), Hi^icb), H2^{b,a). 

In order to have a clearer picture of how these caustics change, let us explain the situation 
when A approaches the singular values c, 6, or a. First, when A — > c~ (respectively, A — > c+), 
the caustic Qx flattens into the region of the coordinate plane — {z — 0} enclosed by 
(respectively, outside) the /oca/ ellipse 

Ql = \{x,y,0) eR':— + -^^l\. (16) 
1^ a — c — c ) 

Second, when A — )■ (resp., A — )■ fe"*"), the caustic Qx flattens into the region of the coordinate 
plane ny — {y — 0} between (resp., outside) the branches of the /oca/ hyperbola 

Ql = l{x,0,z) eR' : ^ I 

I a — — c 

Third, the caustic flattens into the whole coordinate plane tTx = {x = 0} when X ^ a~. 

We recall that not all combinations of caustics can take place. For instance, both caustics 
can not be ellipsoids. The four possible combinations are denoted by EHl, HlHl, EH2, and 
H1H2. Hence, the caustic parameter A = (Ai, A2) belongs to the nonsingular caustic space 

A={ExHi)U {Hi Hi) U{Ex H2) U {Hi x H2), (17) 

where i^i i/i = {A e iJi x iJi : Ai < A2}. For instance, X e E x Hi for trajectories of 
type EHl, which means that Qai is an ellipsoid and Qx^ is a one-sheet hyperboloid. 



The frequency map for billiards inside ellipsoids 



20 



5.2. The extension of the frequency map 

To begin with, we extend the frequency map a; : A — t- to the borders of the four 
components of the caustic space ( 17 ), in the same way that the rotation number was extended 
to the endpoints of the two caustic intervals (10). The extension depends strongly on the 
"piece" of the border under consideration. Hence, we need some notations for such "pieces". 

The set A is the union of three open rectangles and one open isosceles rectangular 
triangle. In total A has eleven edges and eight vertexes. We consider the partitions 

M = A°UA\ A^ = GUi?U5, 

where A^ is the set of edges, A° is the set of vertexes, and S, G and R are the sets formed 
by the four inner edges, the two left edges, and the remaining five edges, respectively. See 
figure |6j We shall see that the frequency map is quite singular (in fact, exponentially sharp) 
at the four edges in S, quite regular at the five edges in R, and it is somehow related to the 
geodesic flow on the ellipsoid Q at the two edges in G. That motivates the notation. 

Next, we shall check that the frequency map of the triaxial ellipsoid Q can be 
continuously extended to the borders of the caustic space in such a way that its values on 
the edges and vertexes can be expressed in terms of exactly six functions of one variable that 
"glue" well. Three of them are the extended rotation functions associated to the three ellipses 
obtained by sectioning Q with the coordinate planes n^, Hy, and tt^. That is, they are the 
functions : [0,b] M, Py : [0, a] — )• M, and : [0, a] — )■ M defined as 

Px(A) = p(A;c,6), Py(A) = p(A;c,a), p^(A) = p(A; 6, a), 

using the notation in ( [T2] ). The other three functions are defined in terms of the former ones 
as follows. Let m = min(A, c) and fn = max(A,c). Let Tx(s) = (A — s){c — s){b — s), 
Ty{s) = (c — s)(A — s){a — s), and Tz(s) = (m — s)(6 — s)(a — s). Then we consider the 
functions : [0, b] — t- M, i/y : [b, a] — )■ M, and : [0, 6] — t- M defined by the identities 



(a 



+2py(A)/ \ 



(6-s)V^ A {s-b)y^) ^PW) 

' 2p,(m) / _^ - , ^\^^ = 0. 



{m- s)^/^J7j Jb {s - m)^/^J7j y^-T^{m) 

Lemma 12. The functions v^, Vy, and have the following properties. 

(i) They are analytic in E U Hi, H2, and E U Hi, respectively. 

(ii) They can be continuously extended to [0, b], [6, a], and [0, b], respectively. 

(Hi) Their asymptotic behavior at the endpoints A G dE U dHi U dH2 = {0, c, b, a} are: 

(a) z/,(A) = 0(Ai/2), asX^ 0+; 

(b) z/,(A) = 0(Ai/2), asX^ 0+; 

(c) z/x(A) = pz(a) + O (1/log |c - A|), as X ^ c; 

(d) z/,(A) = 1/2 + (|A - c|i/2^, asX-> c; 
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(e) z/x(A) = Py{a) + 0(6 — \), as \ ^ h ; 

(f) z/,(A) = p,(c) + O ((6 - A)i/2), asX^ h'; 

(g) z/y(A) = p,(c) + O ((A - 6)^/2), a5 A ^ 6+; and 

(h) z/y(A) = Px(6) + 0(a — A), a5 A — )■ a~. 

Proof. We know that the function Px(A) = p(A; c, 6) is analytic in A, c, and 6, as long as 
< c < 6 and X e E U Hi. Besides, the integrand (a — s)"^(Tx(s))^^/^ is analytic with 
respect to the variable of integration s in the intervals of integration (0,m) and (m, b), and 
with respect to the parameters A, c, b, and a, as long asO<c<6<a and A G EL) Hi. Hence, 
the function z/x(A) = i^x(A; c, 6, a) is analytic in its four variables, as long asO<c<6<a 
and X G E U Hi. The analyticity of Vy and z/^ follows from similar arguments. 

The study of the asymptotic behavior of the functions z/x, Uy, and has been deferred 



to 



Appendix A.9[ Appendix A.10[ and Appendix A.l 1[ respectively. 



□ 



Remark 9. We have numerically observed that and are increasing in E and decreasing 
in Hi, whereas Uy is increasing in H2, but we have not been able to prove it. 

Theorem 13. The frequency map : A — )■ has the following properties. 

( i) It is analytic in A. 

(ii) It can be continuously extended to the border dA, and the extension has the form 



co{0,X2) 


= (0,0) 


for c < X2 


<b, 


w(Ai,6) 


= (Py(^l),Py(^l)) 


for < Ai 


<b, 


w(c, A2) 


= (l/2,p,(A2)) 


for c < X2 


< a, 


w(Ai,a) 


= (Px(Ai),z/x(Ai)) 


for < Ai 


<b, 


u{b, A2) 


= (z/y(A2),Py(A2)) 


for b < X2 


< a, 


w(Ai,c) 


= ('^z(Ai),Pz(Ai)) 


for < Ai 


<c. 


w(Ai, Ai) 


= ('^z(Ai),Pz(c)) 


for c < Xi 


< b. 



(Hi) Its asymptotic behavior at the eleven edges in = G U S U R is: 

(a) uiXi, A2) = K^iX2)Xl^^ + 0(Af as Xi ^ 0+; 



(b) uj{Xi,X2) -a;(c, A2) 

(c) u{Xi,X2) - w(Ai,6) 



K 



,5, 



c, A2)/log |c — All, as Ai — j- c; 



fi;'^(Ai, 6)/ log \b — A2I, a* A2 — )■ 6; ani^ 



(Jj u{X) - w(A^) = 0(A - A^), asX^ X^e R; 
for some analytic functions : HiU H2 ^ and : S ^ 
(iv) Its asymptotic behavior at the eight vertexes in A° is: 



(a) 


uj{Xi 


A2) 


= 0(A;/'), as (Ai,A2)^ 


(0- 


^c+); 






(b) 


uj{Xi 


A2) 


= 0(A;/'), as (Ai,A2)^ 


(0- 


^&); 






(c) 


uj{Xi 


A2) 


= 0(A;/'), as (Ai,A2)^ 


(0- 








(d) 


uj{Xi 


A2) 


= (l/2,p,(c)) + 0(l/log 


■ |c 


— Ai , A2 - 


- c), as (Ai, A2) (c, c+ 


); 


(e) 


uj{Xi 


A2) 


= (1/2, 1/2) + 0(1/ log 1 


c — 


Ai|,l/loj 


5 16 - A2I), as (Ai, A2) -> 


(c, 6); 


(f) 


uj{Xi 


A2) 


= (l/2,pz(a)) + 0(l/lo^ 


l\c 


— Ai , a — 


A2), as (Ai, A2) (c, a" 


); 


(8) 


uj{Xi 


A2) 


= (Py(6),Py(6)) + 0(6- 


Ai, 


l/log|6- 


- A2I), as (Ai, A2) -> (6", 


6); anJ 
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CqCi C2 C3 C4 C5 

EHl: •-f • • • • 

Ai c A2 a 

Co Ci C2C3 C4 C5 

HlHl; • • f-f • • 

c A1A2 a 

Co C1C2 C3C4 C5 

EH2: • • • • • • 

Ai c 6 A2 a 

Co Cl C2C3 C4C5 

H1H2: • • f-f f-« 

c Ai A2 a 



Figure 7. The four possible configurations of the ordered sequence < ci < • • • < C5. Thick 
lines denote intervals of integration. Each one of the displayed configurations illustrates some 
collapse: geodesic flow limit (and type EHl), simple regular collapse (and type HlHl), double 
singular collapse (and type EH2), and double regular collapse (and type H1H2). 



(h) w(Ai, A2) = (px(&), Py(a)) + 0(6 - Ai, a - A2), as (Ai, A2) {h ,a ). 



Proof. Once fixed the parameters a>b>c>Oof the ellipsoid and the couple of caustic 
parameters Ai and A2, we set 

{ci, . . . , C5} = {a, b, c} U {Ai, A2}, Co := < Cl < ■ ■ ■ < C5. 

Four configurations are possible; see figure |7] We said in remark [1] that the frequency is 
analytic in ci, . . . , C5 provided that < ci < ■ • ■ < C5. In particular, this implies that the 
frequency is analytic in the caustic parameter provided it belongs to A. 

The frequency map is expressed in terms of six hyperelliptic integrals over the intervals 
(0, Cl), (c2, C3), and (c4, C5) — represented in thick lines in figure|7] See definition[2| We face 
its asymptotic behavior at the border OA = A°U A^, which requires the study of the asymptotic 
behavior of the six hyperelliptic integrals when some interval defined by the ordered sequence 
< Cl < • ■ ■ < C5 collapses to a point. Therefore, there are exactly five simple collapses. 
The collapse of the first interval is called geodesic flow limit: c\ — > 0+, the collapse of the 
second or fourth intervals is called singular. C2 — Ci — > 0+ or C4 — C3 — t- 0+, and the collapse 
of the third or fifth intervals is called regular: C3 — C2 — ?► 0+ or C5 — C4 — )■ 0+. Thus, regular 
collapses imply that the interval of integration of a couple of hyperelliptic integrals collapses 
to a point; whereas singular collapses imply the connection of two consecutive intervals of 
integration. See figure |7] It is immediate to check that this terminology agrees with the 
partition A^ = G U i? U 5, whereas double collapses — that is, two simultaneous simple 
collapses — correspond to the eight vertexes in A". 

The asymptotic behavior of the frequency map at the eleven edges inA^ = GV} RV} S 
is deduced from several results disseminated through Appendix A[ In short, some technical 
lemmas are listed in Appendix A.l[ some notations are introduced in Appendix A.3[ the 
geodesic flow limit is studied in Appendix A.4[ simple regular collapses are analyzed 
in Appendix A. 5 and simple singular collapses are computed in Appendix A. 6 For instance, 
one can trace the definition of the functions v^, Uy, and z/^ to equation (A.4). The reader 
is encouraged to consult the appendix. Here, we just note that the appendix deals with the 
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general high-dimensional setup, since the computations do not become substantially more 
involved when the dimension is increased. 



The computations regarding the vertexes have also been relegated to Appendix A 



although for the sake of brevity we have written out only the computations for two vertexes. 



Vertex A = (c, h) in Appendix A. 8 — which corresponds to the unique double singular 



collapse — , and vertex A = (6, a) in Appendix A.7 — which correspond to the unique double 



regular collapse. The study of the remaining six vertexes does not require additional ideas. 
For instance, the three vertexes related to the geodesic flow limit can be simultaneously 



dealt with simply by using lemma 22 which ensures that the hyperelliptic integrals over 
(co, ci) = (0, Ai) are 0{\]''^) as Ai ^ 0+. 

Finally, we realize that the extended frequency map w : A — is continuous because 



the extensions "glue" well at the eight vertexes; see lemma 12 For instance, let us consider 
the vertex (6, h). We obtain from the three statements of the theorem regarding this vertex that 

= (Py(6),Py(6)) = (Z/y(6),Py(6)) = ( ( 6 ) , ( C) ) , 

which is consistent: Vyih) = ^^{b) = p^{c) = p{c; 6, a) = p{b; c, a) = Py{b). □ 

Definition 5. The continuous extension u : A ^ R"^ is called the (extended) frequency map 
of the ellipsoid Q. 

The origin of the terminology "geodesic flow limit" can be explained as follows. The 
phase space of the geodesic flow on an triaxial ellipsoid Q C was completely described 
by Knorrer [|29ll . Any nonsingular geodesic on Q oscillates between two symmetric curvature 
lines obtained by intersecting Q with some hyperboloid Qx, A G Hi U H2. The rotation 
number of those oscillations is the quotient 

jmin(fe,A) gds 



p«(A) = -J, T^{s) = -s{X - s)(c - s){b -s){a- s) 



see [fT9l §4. 1]. This rotation number p'^(A) can be continuously extended to the closed interval 
[c, a] with p^{b) = 1. On the other hand, the geodesic flow on the ellipsoid Q with caustic 
hues Q n Qx^ can be obtained as a limit of the billiard dynamics inside Q when its first caustic 
Qai approaches Q; that is, when Ai — )■ 0"^, so that (Ai, A2) — G. Therefore, it is natural to 
look for a relation between the function ) : Hi U H2 and the rotation 

number p'^ : HiU H2 ^ R+. 

Lemma 14. p'^ = K^//«f. Thus, ^2(^1, A2)/a;i(Ai, A2) = p'^(A2) + 0(Ai), as Ai 0+. 



Proof. In Appendix A. 4 we will check that is the unique solution of the linear system 



p-G _LCG \ / G \ / L-G 
-"■01 -"-02 1 I '^l 1 _ I -"-00 

Kf, -Kg 1 I 



where i^g = /;^^+^(T^(s))-i/2sMs, Kg = 2{abc\)-^l\ and {02,03,04,05} = {a, 6,0, A} 
with 02 < 03 < 04 < 05. Therefore, since A G HiLiH2, it turns out that 02 = 0, 03 = min(6. A), 
04 = max(6. A), and 05 = a. Finally, k^/k^ = Kg/ Kg = p^ ■ □ 
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Figure 8. The normalized Jacobian J* : A — [0, 1] of the frequency map for eight different 
ellipsoids. 



5.3. On the Jacobian of the frequency map 

We present the numerical experiments about conjecture [T] stated in section [3] We have 
computed the Jacobian of the frequency map 

J : A ^ M, J(A) := det ( ^(A)^ 

for several ellipsoids, in order to check that it never vanishes. Its visualization close to the four 
inner edges labelled with the letter S in figure [6] has a technical difficulty. To understand this 
fact, one can look at the graph of the rotation number p(A) shown in figure [3j The derivative 
p'(A) explodes at A = 6, which would make difficult its visual representation. The problem is 
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worse in the spatial case, because the frequency map has the same kind of "inverse logarithm" 
singularity at the four inner edges instead of at a single point. 

We overcome the visualization problem by representing the normalized Jacobian 

J. : A ^ [0, 1], J.(A) = (1 - exp(-| J(A)|))i/^ 

The exponential function is intended to cancel the exponentially sharp behavior of the 
Jacobian at the inner edges. The exponent 1/4 has been chosen by trial and error to obtain 
more informative plots. The normalized Jacobian ranges over the interval [0, 1]. We note that 
J* = <S=^ J = and J^, = 1 <^ \ J\ = oo. The results are shown in figure [S] In the upper left 
comer, we have displayed the parameter space P introduced in ([15]), and sketched in figure |5] 
We study the eight ellipsoids that correspond to the eight points in P labelled from 1) to 8). 
In particular, we have chosen at least one sample of each "kind" of ellipsoid: 1) standard, 
2) almost spheric, 3) standard, 4) almost prolate, 5) almost oblate, 6) close to a segment, 7) 
close to a flat solid ellipse, and 8) close to a flat circle. The color palette is a classical one: 
cold colors represent low values, hot colors represent high values. The neighborhood of the 
inner edges is always a "hot" region; that is, the Jacobian is always big on that region. On the 
contrary, the Jacobian tends to zero close to the hypotenuse of the Hi ® Hi region. This can 
be seen from a symmetry reasoning. Furthermore, the Jacobian never vanishes, not even in 
the cases 7) and 8), which correspond to almost flat ellipsoids. 



5.4. The range of the frequency map 



We recall that if the two conjectures stated in subsection 3.2 hold, then the components of the 
frequency map are ordered as stated in (|9]). Thus, the range of the frequency map should be a 
subset of the frequency space 

Vt = {iuJi,U2) G : < ^2 < wi < 1/2} . 

We visualize in figure |9] how each edge of the caustic space is mapped onto the frequency 
space. All the depicted curves have been numerically computed from exact formulae given 



in theorem 13 We have represented the caustic space A at the left side, and the frequency 
space ri at the right side. Each colored segment in the caustic space is mapped onto the curve 
of the same color in the frequency space. The black segment in A — which represents the 
geodesic flow limit — is mapped onto the origin O = (0, 0). The point (c, b) is mapped onto 
A = (1/2, 1/2). The images of the magenta and blue segments are folded at this point A. 
Henceforth, [AB] stands for the segment with endpoints A and B, and A[ABC] stands for 
the interior of the triangle with vertexes A, B, C . We see that uj{E x Hi) is enclosed by 
the magenta segment [0A\, the blue segment and a red smooth curve from Bi to O; 

uj{Hi ® Hi) = A[ABiCi]; u{E x H2) is enclosed by the magenta segment [OA], the blue 
segment [Ai?2], and a cyan smooth curve from B2 to O; and co{Hi x H2) is enclosed by the 
magenta segment [CiA], the blue segment [v4i?2], a brown smooth curve from i?2 to D, and a 
green smooth curve from D to Ci. 




The points A, Bi, B2, Ci, and D can be explicitly expressed in terms of the parameters 
of the ellipsoid Q. Let < g^, Qy, Qz, Q* < 1/2 be the quantities defined by 



sin^ Tig^ = c/b, sin^ HQy 



= p{c;b,a), (18) 



we 



c/a, sin ng^ = b/a 

where p(A; b, a) is the rotation number ( 12 1. From the formulae contained in theorem 13 
get that A = (1/2, 1/2), Bi = (1/2, g,), B2 = (1/2, g,), Ci = {g,, g,), and D = (&, g^). 
We note that D E ^l, since gy < g^. In fact, gy < g^ and gy < g^, < g^, although we do not 
have a rigorous proof of the inequalities involving g^:. The four quantities defined in ( 18 1 can 
be interpreted in terms of the restriction of the billiard dynamics to suitable planar sections of 
the original ellipsoid. For instance, g* is the rotation number of the trajectories contained in 
the section by the plane tt^ whose caustic is the focal ellipse ([16]). 

We give now some numerical estimates on the size and the shape of the four ranges. 



Numerical Result 2. Let < gx, Qy, Qz, (>* < 1/2 be the quantities defined in {18). Let 
g = max(f?„&). Let O = (0,0), A = (1/2,1/2), B^ = (1/2, f?,), B^ = (1/2, f?,), 
B = (1/2, g), Ci = {g^, g^), C2 = (&, &), C = {g, g), and D = {g^,, gy). Then: 

(i) A[ABjCj] C uj{E X Hj) C A[ABjO]forj = 1,2; 

(ii) uj{Hi ® Hi) = A[ABiCi]; and 

(Hi) A[ABC] C uj{Hi X H2) C uj{E x H2) n {{ui,U2) e Q : ui > g^, UJ2 > gy}. 

Next, we enlighten some practical consequences of these estimates. To begin with, let us 
present four simple criteria to decide if the ellipsoid has billiard trajectories of frequency 
a;° = {uJi,u}^) e Q and of caustic type EHl, HlHl, EH2, or H1H2. Compare with 
the criterion for the existence of billiard trajectories inside an ellipse with rotation number 
p° G (0, 1/2) and a caustic hyperbola given in ( 13 1. 
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Proposition 15. If numerical result^holds, then the following criteria can be applied. 

(i) Ifu^ > p{c; b, a), then G uj{E x Hi). Ifuj^/2u'l < p{c; b, a), then a;° ^ u{E x Hi). 

(ii) oj^ G uj{Hi (g) Hi) if and only z/a;^ > p(c; b, a). 

(Hi) Ifb < asin^vru;"' ^ ^(-^ ^ H2). Ifb > a sin^(7rw2/2'^i). then w° ^ cj(i? x H2). 

(iv) If > p{c;b,a) and b < asin^TTW^, then G uj{Hi x if2)- //''^i < p{c;b,a), or 
c> a sin^ vrw^, orb>a sin^(7ra;2/2a;5'), f/ien o;*^ w(-f^i x /72)- 

Hence, there exist billiard trajectories of the four caustic types when is big enough: 
IJJ2 > p{c;b,a) and sin^ ttw^ > b/a. On the contrary, there does not exist any of such 
trajectories when w^/w? is small enough: ^2/1^1 < 2p(c; b, a) and sm'^{Tcu2/2ui) < b/a. 

Proof The first and third criteria follow from A[ABjCj] C uj{E x Hj) C A[ABjO], 
j = 1,2. The second one follows from the identity uj{Hi (g) Hi) = A[ABiCi]. The last 
one follows from the the last item of numerical result [21 □ 
We can also understand how the range of the frequency map depends on the shape of 
the ellipsoid. It suffices to see how the quantities g^, Qy, Qz, and depend on the parameters 
< c < 6 < a. On the one hand, if the ellipsoid flattens — that is, if c decreases, but a and b 
remain fixed — , then f?* decreases, so to (E x Hi) anduj{Hi®Hi) expand. Indeed, both ranges 
tend to cover the whole space for flat ellipsoids: c — > 0"^, whereas they collapse to the empty 
set for prolate ellipsoids: c — )■ 6 . On the other hand, if the ellipsoid becomes more oblate 
— that is, if b increases, but a and c remain fixed — , then increases, so uj(E x H2) contracts. 
Indeed, uj{E x H2) tends to cover f2 for "segments": b — )■ 0+, but collapses to the empty set 
for oblate ellipsoids: b . The behavior of to (Hi x H2) is more complicated, because its 



vertex D = (^x, Qy) can be at any point of the frequency space Vl; see (18). Anyway, if the 
ellipsoid becomes spheric — that is, c and b approach a — , then and tend to one half, so 
Bj tends to A and the four ranges collapse to the empty set. This means that the more spheric 
is an ellipsoid, the poorer are its four types of nonsingular billiard dynamics. Some of the 



criteria stated in propositions 15 and 17 quantify this general principle. 



The ranges of the frequency map for eight different ellipsoids are shown in figure 10 



In the left upper picture, we have again marked the ellipsoids as points in the parameter 



space ( 15 1. The image sets u{E x Hi), uj{Hi®Hi), lu{E x H2), and uj{Hi x H2) are depicted 
in yellow, green, magenta, and blue, respectively. The transparency allows to visualize 
simultaneously all four sets. We can check all their properties stated in numerical result |2} 
together with the ones regarding their dependence on the shape of the ellipsoids. Blue dots 
correspond to rational frequencies with small common denominators. 

In the previous paragraphs the range of the frequency map has been described by mixing 
analytic formulae and numerical computations, but some properties can be justified. The 
following proposition is an example. 

Proposition 16. If conjecture\ljholds, then uj{Hi ® Hi) = A[ABiCi] C uj[E x Hi), and 
uj : Hi ® Hi ^ A[ABiCi\ is a global dijfeomorphism. Besides, G ijj{Hi Hi) if and 
only ifuj2 > g* '■= p{c; b, a). Finally, limc„i.o+ = Q and limc^fe- A[ABiCi] = 0. 
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Figure 10. Ranges of the frequency map for eight different ellipsoids. 



Proof. If f/ = Hi Hi, then X = dU is the triangle with vertexes A = (c, b), Bi = (c, c), 
Ci = {b,b). Using the formulae for the extended frequency map established in theorem 1 3 , we 



get that = [ABi], uj{[BiCi]) = [BiCi], and uj{[CiA]) = [CiA]. Thus, Y = io{X) 

is the triangle with vertexes A, Bi, Ci. In particular, X and Y are Jordan curves, so the 



frequency map a; : f/ — t- M verifies the hypotheses of lemma 27 in Appendix B Hence, 



uj{U) = A[ABiCi] and u : U ^ ^{U) is a global diffeomorphism. 

In order to prove the strict inclusion A[ABiCi] C uj{E x Hi), it suffices to see that the 
red curve from O to Bi in the right picture of figure |9] is strictly below the yellow segment 
[-BiCi]. And this is equivalent to prove the inequality 

Pz(Ai) < Pz(c), VAiG(0,c), 
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Type 


nil 


1712 


UJl 


UJ2 


EHl 


Crossings of tTz 


Half-turns around 


z-oscillation 


z-rotation 


EH2 


Half-turns around 


Crossings of tTx 


x-rotation 


x-oscillation 


HlHl 


Touches of Qx. 


Half-turns around Oz 


(Hl-oscillation)/2 


z-rotation 


H1H2 


Crossings of vTy 


Crossings of tTx 


y-oscillation 


x-oscillation 



Table 1. Geometric meaning of the winding numbers and the frequency vector when Q C M'^. 



due to the formulae for the extended frequency map contained in theorem 13 This inequality 



was proved in proposition 11 Finally, we note that limc_>.o+ g* = and limc^^- = 1/2; 



see the second item of proposition [10} □ 



5.5. Geometric meaning of the frequency map 

Let mo, rui , m2 be the winding numbers of a periodic billiard trajectory of type EHl . Then mo 
is the period. Besides, according to remark|2} mi and 1122/2 are the number of times along one 
period that the trajectory crosses the coordinate plane = {z = 0} and the number of times 
along one period that it rotates around the coordinate axis = {x = y = 0}, respectively. 
Therefore, the components of the frequency map have the following geometric meaning: 
coi = mi/2mo is the number of oscillations around tt^ per period, whereas 102 = 1112 /2mo 
is the number of rotations around per period. Thus, it is quite natural to say that coi is the 
z-oscillation number and U2 is the z-rotation number of the trajectory. 

As in the planar case, these interpretations are extended to quasiperiodic trajectories. If 
X = {Xi, X2) E E X Hi, then Qai is an ellipsoid, is a one-sheet hyperboloid, and 

io{X) = lim {nk,lk)/k, 

fc— >+oo 

where Uk (respectively, Ik) is the number of oscillations around vr^ (respectively, number of 
rotations around a^) of the first k segments of any given trajectory with caustics Qai and Qx^ . 

The billiard trajectories of other types can be analyzed following similar arguments. The 
results are listed in table [T] and can be checked by visual inspection; see figure 13 



Finally, we stress a point already commented in remark [3| If the trajectory is of type 
HlHl — that is, if both caustics are one-sheet hyperboloids — , then the winding number mi 
is the number of (alternate) tangential touches with the caustics, so ui = mi/2mo is half 
the number of oscillations between the one-sheet hyperboloids per period. In that situation, 
we call 2wi the HI -oscillation number of the trajectory. In particular, it can happen that 
mow ^ Z^. For instance, if the winding numbers are mo = 4, mi = 3, and m2 = 2, the 
period is four, but a; = (3/8, 1/4). 

5. 6. Bifurcations in parameter space 

We want to determine the ellipsoids that have billiard trajectories with a prescribed frequency 
and with a prescribed caustic type. We recall that each ellipsoid is represented by a point in 
p = {(6, c) G ; < c < 6 < 1}, because a = 1. Let P°, P°. ^3 ^ and P° be the four 
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regions of P that correspond to ellipsoids with billiard trajectories of frequency uj^ and caustic 
type EHl, HlHl, EH2, and H1H2, respectively. Their shapes are described below. 

Numerical Result 3. Once fixed any frequency vector uj^ = (w°;'^2) ^ let = 6° = 1- 
63 = 64 = sin^(7rci;2/2ci;5'), c1 = = (3^ = sin^ ttcj^/ si'^^ '^^v <^ndc2 = c\ = sin^ ttw^. Then 

P^ = {{h,c)eP:h<hl c<g%h)], 1<J<4, 



for some continuous functions (7° : [0, iPA — )■ M 5Mc/z 



.0. 

1' 



.0. 

-2' 



(«j is concave increasing in [0, 1], < Qiip) < hfor all b G (0, 1), and (??(!) = c' 
(ii) §2 is concave increasing in [0, 1], < g2ip) < Qiip) for all b E (0, 1), and (72(1) 
(Hi) is the identity in [0, C3], concave decreasing in [c^, 63], and g^ib^) = 0; and 
(iv) g^ is increasing in [0,(3^], concave decreasing in [f3^,b1], c^b/f3^ < g^{b) < b for all 
b G (0,/30), < glib) < glib) for all be (/J", 6°), gliPl) = 4 and gl{bl) = 0. 

Remark 10. Numerical result [T] follows from numerical result [3] just by choosing suitable 
rational frequency vectors: tu^ = (2/5, 1/5) in the cases EHl and EH2, w° = (3/8, 1/4) in 
the case HlHl, and w° = (1/3, 1/6) in the case H1H2. We stress that inequality gl < g^ in 
numerical result [T] and inequality g^ < g\ in numerical result |3] are not contradictory, because 
the first one refers to two different frequency vectors: (2/5, 1/5) and (3/8, 1/4). 

Remark 1 1. We have numerically checked that g^ is not concave in [0, (3l]. 

Remark 12. Inclusions P2 C P° and P4 C P3 — and so, inequalities (^) < 9i{^) 
glib) < glib) — are in direct agreement with inclusions a;(iJi ifi) C coiE x Hi) and 
uiHi X H2) C uiE X H2) mentioned in numerical result[2j 

Some bifurcations curves corresponding to the graphs of the functions g^ : [0, — M 
are presented in figure 1 1 On top of this figure we consider the eight rational frequencies with 
the smallest denominators. The inclusions P2 C P° and P4 C P3 can be easily visualized, 
since all dashed curves are below their continuous pairs. On the bottom, we depict the 



bifurcation curves associated to the rational frequencies marked with blue dots in figure 10 



We have needed a multiple precision arithmetic to compute the bifurcation curves close to 
some of their endpoints, since the involved root-finding problems become quite singular at 
them. The programs have been written using the PARI system [5]. 

Next, we describe four more criteria to decide if an ellipsoid has billiard trajectories of a 
given frequency. They are similar to the four ones established in proposition [T5| 

Proposition 17. If numerical result^holds, then the following criteria can be applied. 

(i) Ifc < c\b, then G ujiE x Hi). If c > c^a, then uiE x Hi). 

(ii) Ifc < c^b, then G w(i/i (g) Hi). Ifc > c^a, then ^ w(/7i (g) Hi). 

(Hi) 7/(6° - c°)c < c°(6°a - b), then co^ G uiE x H2). Ifb > b^a, then co^ ^ uiE x H2). 

(iv) If(3lc < min (c^6, /3°6°a + (c° - 6°)6), then w° G uiHi x H2). Ifb> bla or c > c\a, 
then w° ^ i^iHi x H2). 
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2/5, 1/5 


1/3, 1/6 


2/7, 1/7 


3/7, 1/7 


3/7, 2/7 


1/4, 1/8 


3/8, 1/8 


3/8, 1/4 




Figure 11. Some bifurcation curves in the parameter space P. 



Proof. From numerical result |3 

V\ = Vl = (1,0), ro = (60,0) 



and A° 



we get that := /\[OT]A]] c P°, where O = (0,0), 
, ro = (60,0), A? = (l,c?), AO = (l,cO), AO = (cO,cO), 
(/30,cO). It is straightforward to check that a point (6, c) G P belongs to the 
triangles T^, T^, T^, and if and only if c < c?6, c < c^6, (6^] - c|])c < c!](6|] - b), and 
/^Oc < min {c^b, l3^b^ + (cO — 64)6), respectively. This proves the first part of each criterion 
for a = 1. To prove the general case, it suffices to take into account that its formulae are 
homogeneous in the parameters a, b, c. 

The second parts follow from similar arguments. For instance, Pj* fl {c > cO} = for 
j = 1, 2, since gj{b) are increasing in [0, 1] and (7° (1) = c^- □ 
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Remark 13. Proposition 15 has been obtained by fixing the ellipsoid and looking at the 
frequency space. On the contrary, proposition 17 has been derived by fixing the frequency 
vector and looking at the parameter space. Of course, both approaches are equivalent, but 
their criteria are slightly different. The second ones are computationally simpler, because 
they do not involve any elliptic integral. 

Although the description of the regions has a strong numerical component, some 
results can be proved. The following proposition is an example. 

Proposition 18. Let = (cuJ'jCJ^) G be a fixed frequency vector If conjecture^holds, 
then C P°, P2 only depends on u^, and P^ = {{b,c) e P : < b < 1, c < g^{b)} for 
some increasing analytic function g2 : (0, 1) — )■ M such that < g2{b) < bfor all b G (0, 1). 
Besides, limfe^i- g2{b) = sin^ vrcj^- Finally, 

36/(1 + 6 + 2^1-6 + 62) foru;^ = l/3, 

6/(1 + 6) foru;^ = l/4. 

Proof If conjecture [T] holds, then uj{Hi ® Hi) = A{ABiCi) C u{E x Hi); see 
proposition 16 Therefore, because 

(6, c) G P^ ^ w° G uj{Hi (8) Hi) C uj{E x Hi) (6, c) G Pf. 

On the other hand, since A = (1/2, 1/2), Bi = (1/2, f?*) and Ci = (f)*, g^:), we deduce that 

(6, c) G P^ ^ w° G uj{Hi ® Hi) = A{ABiCi) ^ tu^ > := p{c; 6, a). 

We know that the rotation function p(A) = p(A; b,a) is increasing in (0,6), p(0) = 0, and 
p(6) = 1/2. Hence, the function g2 : (0, 1) M, < gf^l^) < b, is implicitly defined by 

p((?°(6);6,a)=c^°, (19) 

where G (0,1/2) and a = 1 are fixed parameters. Analyticity of g2 follows from 
the Implicit Function Theorem, since conjecture [T] also implies that (9ip(A; 6, a) 7^ for 
all A G (0,6) U (6, a). Indeed, this derivative is positive in (0,6) and negative in (6, a), 
because p(-;6, a) is increasing in (0,6) and decreasing in (6, a). Besides, we know that 
(92p(A;6, a) = dip{b;\,a) from the symmetry p{X]b,a) = p(b]\,a); see (12). Hence, by 
differentiating equation ( 19 1 with respect to 6 and setting c = (7° (^) ^ (O5 we get that 

(92)' (b) = -d2pic;b,a)/dip{c;b,a) = -dip{b; c, a) /dip{c;b, a) > 0. 

Using proposition[To[ we know that lim;,_j.a- sin^ 7rp(c; 6, a) = hm^.^^- sin^ 7i"p(6; c, a) = c/a. 
Thus, we deduce lim^^i- (^^l^) = sin^ vrw^' since a = 1. 

Finally, we must find the values c G (0, 6) such that p(c; 6, a) is equal to 1/3 or 1/4. That 
is, we must find the values of c G (0, 6) such that the billiard trajectories inside the ellipse 
{x2/a + ?/2/5= 1} with caustic {x2/(a — c) +?/2/(5 — c) = 1} have period three or four. This 
is an old result that goes back to Cay ley dH. For instance, c = a6/(a + 6) in the four-periodic 
case. The value for the three -periodic case was given in equation ([14]). □ 



The fact that P2 only depends on UJ2 can be visualized on top left in figure 1 1 The two 
dashed curves with = 1/8 coincide, as well as the two ones with = 1/7. 




5.7. On the ubiquity of almost singular trajectories 



In figure 12 we have superposed the edges and borders (drawn in light colors) already 
displayed in figure |9| and some new segments and curves (drawn in heavy colors). In the 
caustic space, these new segments are close to the original edges. To be precise, the distance 
between them and the edges is equal to c/100 = 4.6 ■ 10^'^. Nevertheless, the images of the 
black, magenta, and blue ones are far from their corresponding borders in the frequency space. 
This phenomenon seems stronger on the magenta and blue borders. It has to do with the fact 



that, as stated in theorem 13 the frequency map has an inverse logarithm singularity at the 
blue and magenta edges of the caustic space. Therefore, one must be exponentially close to 
them, just to be close to their images. On the other hand, the frequency map has a squared 
root singularity at the black edges of the caustic space. Thus, one must be quadratically close 
to them, just to be close to the origin in the frequency space. 

We deduce from this phenomenon that billiard trajectories with some almost singular 



caustic are ubiquitous. Let us describe a quantitative sample of this principle using figure 12 
Let T be the triangle delimited by the yellow, blue and magenta thin segments that are close 
to the edges of Hi ® Hi. It turns out that the area of uj{Hi ® Hi) is approximately 16 
times the area of uj{T). Hence, if we look for billiard trajectories of type HlHl inside Q 
with a random frequency in uj{Hi ® Hi), their caustic parameter A = (Ai, A2) shall verify 
min(| Ai — c|, | A2 — 6|) < 4.6 • 10"^ with a probability approximately equal to 94%. It suffices 
to note that 15/16 = 0.9375. 
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(a,5,c) = (1, 0.4, 0.2), mo = 5 (a, 6, c) = (1, 0.4, 0.2), toq = 6 
(Ai,A2) = (0.1995,0.7630), (wi,W2) = (2/5,1/5) (Ai,A2) - (0.2191, 0.9234), (wi, cja) = (1/3,1/6) 





(Ai,A2 



(,a,6,c) = (1, 0.49, 0.25), mo = 5 
= (0.2316,0.2603), (^1,^2) = (2/5, 1/5) 



(a. 6, c) 
(Ai,A2) 



= (1, 0.81, 0.36), 
= (0.3696,0.6087), 



Too = 4 

(wi,W2) = (3/8, 



1/4) 



Figure 13. Examples of symmetric nonsingular billiard trajectories with minimal periods. 
Lines in red represent the particle's trajectory. Lines in green and yellow correspond to 
the intersections of the original ellipsoid with the caustic 1 -sheet and 2-sheet hyperboloids, 
respectively. In the cases EHl and EH2, the caustic ellipsoid is also depicted. 



5.8. Examples of periodic trajectories with minimal periods 

We have numerically computed some symmetric periodic trajectories to check that the lower 
bounds stated in theorem [T] are optimal; see figure [13] All these trajectories are almost 
singular. Concretely, c — Ai ~ 5 • 10^"^ in the case EH2; \i — c c:^ 2 ■ 10~^ in the case 
H1H2; c — Ai ~ 2 ■ 10~^ and A2 — c ~ 10~^ in the case EHl; and Ai — c ~ 10"^ in the case 
HlHl. Of course, we did not look for almost singular trajectories, but we got them anyway. 

Considering the values given in figure [13} and bearing in mind table [Tj we have 
(mi, 7712) = (4,2) for the EH2 trajectory, so it performs two turns around the coordinate 
axis Ox and crosses twice the coordinate plane n^. As well, (mi, 7712) = (4,2) for the 
EHl trajectory, meaning four crossings with tt^ and just one turn around a^. Again, we 
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have {mi, 1712) = (4,2) for the H1H2 trajectory, meaning four crossings with TTy and two 
crossings with n^. Finally, {mi, 1112) = (3, 2) for the HlHl trajectory, which corresponds to 
three tangential touches with each of the caustics and a single turn around Oz. Each of those 
geometric interpretations has been verified on the corresponding trajectory. 

6. BUliard inside a nondegenerate ellipsoid of ]R"+^ 

We describe briefly the high-dimensional version of some of the analytical results already 
shown in the spatial case. We denote again the nondegenerate ellipsoid as in ([T]) and the 
nonsingular caustic space as in ([2]). 

By analogy with the spatial case, we consider three disjoint partitions: 

dA = U^-(JA^ A''-^ = GURUS, S = U]^^Sj. 

With regard to the first one, A'^ is the /c-dimensional border of A. That is, A° is the set of 
vertexes, A^ is the set of edges, A^ is the set of faces, and so on. The second one mimics 
the distinction among geodesic flow limits, simple regular collapses, and simple singular 
collapses already seen in the previous section. For instance, G = {\ E A"^^ : Ai = 0}. The 
asymptotic behavior of the frequency map in each one of these three situations is expected to 
be dramatically different; see the theorem below. The last partition labels the component of 
the caustic parameter that becomes singular: Sj = {\ G A"~^ : = aj}. Besides, given any 
caustic parameter A G A we shall denote by A'^^ G Sj the caustic parameter obtained from A 
by substituting its j-th component with aj. Finally, we introduce the {n — 1) -dimensional set 

a = {(A2, . . . , A„) G : (0, A2, . . . , A„) G G}, 

which turns out to be the nonsingular caustic space for the geodesic flow on the ellipsoid. We 
note that Si = {c} x {Hi U H2), S2 = {EU Hi) x {b}, and G^ = HiU H2 with the notations 
used in the previous section for triaxial ellipsoids of M^. 

Theorem 19. The frequency map a; : A — t- has the following properties. 

( i) It is analytic in A. 

(ii) It can be continuously extended to the border dA, the extended map being as follows: 

(a) It vanishes at G; 

(b) One of its components can be explicitly written as a function of the rest at R; 

(c) Its first component is equal to 1/2 at Si; 

(d) Its l-th component is equal to the {I — l)-th component at Si for 2 < I < n; and 

(e) Its "free" components are an {n — 1) -dimensional frequency of the billiard inside 
the section of the original ellipsoid by a suitable coordinate hyperplane at RU S. 

Besides, the restriction of the continuous extended map to any of the k-dimensional 
connected components of A^, 1 < k < n — 1, is analytic. 

(Hi) Its asymptotic behavior at A"^^ = G U S U R is: 

(a) uj{X) = K^{X2, Xn)Xy^ + 0(Af as Ai ^ 0+; 

(b) u{X) - u{X^^) X K^{X^^)/\og \aj - Xj\, as Xj 
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(c) u{X) - uj{\^) = 0(A - A^), asX^X^e R; 
for some analytic functions : ^ M" and k,^ : S ^ M". 

Proof. It follows from the same arguments and computations that in the spatial case. The 
arguments are not repeated. The computations with hyperelliptic integrals have been relegated 
to [Appendix A[ □ 
We recall that, once fixed the parameters ai, . . . , a„+i of the ellipsoid and the caustic 
parameters Ai, . . . , A„, we write the 2n + 1 positive numbers 

{Ci, . . . , C2n+l} = {Ol, • • • , Cln+l} U {Ai, . . . , A„} 

in an ordered way: cq := < ci < ■ ■ ■ < C2n+i- Then the frequency uj{\) is defined in terms 
of some hyperelliptic integrals over the intervals {c2j, C2j+i). If two consecutive elements of 
{co, . . . , C2n+i} collide, then cu{X) is, a priori, not well-defined. Thus, it is natural to ask: How 
does a; (A) behave at these collisions? 

In the previous theorem we have solved this question at the set A"~^ = GURUS, which 
covers just the geodesic flow limit: ci — 0"*", the n simple regular collapses: C21+1, C21 — )■ c* 
for some /, and the n simple singular collapses: 021-1,021 0* for some /. But there are many 
more (multiple) collapses, from double ones to total ones. Double collapses correspond to the 
set A"~^. Total collapses have multiplicity n, so they correspond to set of vertexes A°. 

We believe that it does not make sense to describe the asymptotic behavior of the 
frequency map at all of them, since the behavior in each case must be the expected one. 
In order to convince the reader of the validity of this claim, we end the paper with a couple of 
extreme cases. 

As a first example, let us consider the vertex A = (ai, . . . , a„) G A°. It represents 
the unique total singular collapse, because it is the unique common vertex of the 2" open 
connected components of the caustic space: 

n 

n A. = {A} = n^r 

o-e{o,i}" i=i 



It IS 



Using that the point A belongs to all the closures Sj, from theorem 19 we get that uj{X) 
(1/2, ... , 1/2). Which is the asymptotic behavior of oj at this vertex? In Appendix A.8 
proved that 

w(A) = u{\) + 0(1/ log \ai - All, . . . , 1/log |a„ - A„|), A -> A. 

This behavior is singular in the n caustic coordinates, as expected. 

On the contrary, the vertex A = (02, ... , a„+i) G A° represents the unique total regular 
collapse, so we predict a regular behavior in the n caustic coordinates. In [Appendix A.7|we 



show that a; (A) = (wi, . . . ,Un), where the limit frequencies < uj < 1/2 are defined as 
sin^ TTUJj = ai/ttj+i, and the asymptotic behavior is 

a;(A) =a;(A) + 0(A-A), A ^ A". 

Once more, the frequency map has the expected behavior. 
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7. Conclusion and further questions 

We studied periodic trajectories of billiards inside nondegenerate ellipsoids of ]R"+^. First, we 
trivially extended the definition of the frequency map u to any dimension, we presented two 
conjectures about uj based on numerical computations, and we deduced from the second one 
some lower bounds on the periods. Next, we proved that uj can be continuously extended to 
any singular value of the caustic parameters, although it is exponentially sharp at the "inner" 
singular caustic parameters. Finally, we focused on ellipses and triaxial ellipsoids, where we 
found examples of trajectories whose periods coincide with the previous lower bounds. We 
also computed several bifurcation curves. Despite these results, many unsolved questions 
remain. We indicate just four. 

The most obvious challenge is to tackle any of the conjectures, although it does not look 
easy. We have already devoted some efforts without success. We believe that the proof of 
any of these conjectures requires either a deep use of algebraic geometry or to rewrite the 
frequency map as the gradient of a "Hamiltonian"; see [|45l §4]. 

Another interesting question is to describe completely the phase space of billiards inside 
ellipsoids in M""*"^ for n > 2. A rich hierarchy of invariant objects appears in these billiards: 
Liouville maximal tori, low-dimensional tori, normally hyperbolic manifolds whose stable 
and unstable manifolds are doubled, et cetera. For instance, the stable and unstable invariant 
manifolds of the two-periodic hyperbolic trajectory corresponding to an oscillation along the 
major axis of the ellipsoid were fully described in lfT4ll . 

Third, we plan to give a complete classification of the symmetric periodic trajectories 
inside generic ellipsoids [lOJ. To present the problem, let us consider the symmetric periodic 
trajectories inside an ellipse displayed in figure |4j On the one hand, the three-periodic 
trajectory drawn in a continuous red line has an impact point on (and is symmetric with respect 
to) the X-axis. On the other hand, the four-periodic trajectory drawn in a dashed green line 
has a couple of segments passing through (and is symmetric with respect to) the origin. It is 
immediate to realize that there do not exist neither a trajectory with a hyperbola as caustic like 
the first one, neither a trajectory with an ellipse as caustic like the second one. The problem 
consists of describing all possible kinds of symmetric periodic trajectories once fixed the type 
of the n caustics for ellipsoids in M""*"^. Once these trajectories were well understood, we 
could study their persistence under small symmetric perturbations of the ellipsoid, and the 
break-up of the Liouville tori on which they live. Similar results have already been found 
in other billiard frameworks: homoclinic trajectories inside ellipsoids of ]R"+^ with a unique 
major axis [[8l, and periodic trajectories inside circumferences of the plane II361I . 

Finally, we look for simple formulae to express the caustic parameters Ai, . . . , A„ that 
give rise to periodic trajectories of small periods in terms of the parameters ai, . . . , a„+i of 
the ellipsoid. As a by-product of those formulae, one can find algebraic expressions for the 
functions g*j{h) that appear in numerical resultflj This is a work in progress ll37l . 
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Appendix A. Computations with hyperelliptic integrals 

Appendix A. 1. Technical lemmas 

Lemma 20. Let fe e C°([q;, ^]) be a family of Junctions such that fe — /o + 0(e) in the 
C^-topology. Then 



J a \/{s - a)iB - S) J a 



+ 0(e). 



Proof \Ie-Io\ < |/e-/o|c"([a,/3]) Xf((s-Q;)(/5-s))-^/Ms = 7r|/e-/o|c«([a,/9]) = 0(e)- n 
Lemma 21. Let f e C^([m, M]) with m < a < P < M and e = P - a. Then 

r , = = nf{a) + 0(e) = 7r/(/3) + 0(e), e ^ 0+. 

Proof. Using the Mean Value Theorem for integrals, we get that there exists some Sq e [a, (5] 
such that the integral is equal to /(sq) J^{{s — a){/3 — s))~^/^ds = 7r/(so). □ 

Lemma 22. Let / e C^([0, M]) with < e < M. Then 

I, = /' 4^ = 2/(0)e^/2 ^ Q^^3/2)^ ^ ^ Q+_ 
Jo Ve-s 

Proof /,= [-2(e-s)i/2/(s)];:; + 2/;(e-s)VY(s)ds = 2/(0)eV2 + O(e3/2). □ 

Lemma 23. Let f e C\[a,l3]). Setr] = f{a) log(4/3-4«) + Jj^ {s - a)-\f (s) - f {a))ds, 
e = (s - a)-3/2(/(.) - f(a))ds, /X = /(/3) log(4/3 - 4a) + £^(/3 - .)-H/(^) - /(/3))ds, 
anJt^ = (/3 - s)-3/2(/(s) - /(/3))ds. T/ien 

/■'^ ffs)ds 

/ ^==44;^=== = -/(«) log e + 77 + 0(e log e), e ^ 0+, 
Ja \/{s + e — a)[s — a) 

f , ^ = 7r/(a)e-V^ + i + 0{e^'\ e ^ 0+, 

r , ^^'"^"^^ =-/(^)loge + At + 0(eloge), e ^ 0+, 

L V(/3 + e-s)(/3-s) -^^^^ 6 ^ ^ 6 ^' 

./a (/3 + e-s)V/:^-s 
r/ze yzri'? ( respectively, last) two estimates also hold when f has a singularity at s = {3 
(respectively, at s = a), provided f G L^([a, /?]). 
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Proof. We split the first integral as /g = 77 + /e — /g, where fj — J^{s — a)~^{f{s) — f{a))ds 
is a constant, and 

By performing the change = s — a in the integral /g, we get that 

Jo Vx^ + e L ^ yix=o 

where fj = f{a) log(4/3 — Aa) is another constant. Thus, to get the first formula with constant 
rj — fj + fj it suffices to see that 1^ — 0(e log e). 

Once fixed some 7 e (a, we decompose the integral le as the sum + Xg, where 
J. = r/(s)r,(s)ds, k, = /''/(s)r.(s)ds, and 



First, we consider the interval [a, 7]. Then |/|oo = niax{|/(s)| : a < s < 7} < 00 and rg(s) 
is positive in [a, 7]. Set 5 = 7 — 0;. Using again the change x'^ — s — a, we see that 



J a M s + e-a Jo 



x=Vs 

S- ' ... 



\/ x^ + e 

xvx^ + e + elog x + vx^ + e = — lege + 0(e) = O(eloge). 

\ /Ja;=0 2 

Concerning the other interval, we note that rjys) is positive and decreasing in [7, ^]. Hence, 

max{|rg(s)| : 7 < s < ^9} = tJ^^), and so \K^\ < r^i'j) |/(s)|ds = 0(e). This ends the 
proof of the first formula. 

We split the second integral as Lg = ^ + Lg — Lg, where ^ is the constant given in the 
statement of the lemma, and 

Lg= r ^^"'^^^ Lg= r M^jM(i-^^^)ds. 

Ja (s + e - a)\/s - a ^ Ja (s - q;)3/2 \^ s + e-aj 

By performing the change x = s — a in the integral Lg, we get that 

Lg = ['"^ r^^^ = 2f{a)e-'/' [atan (e-^/^xV^)]^:^"" = 7r/(a)e-V2 + 0(e^/^). 
Jo [x + e)^^/x 

Thus, to get the second formula it suffices to see that Lg = 0(e^/^), which follows from 
similar computations than the ones above. 

The last formulae are obtained by performing the change of variables s — a — P — tin 
the former ones. □ 

Corollary 24. Let f e C^{[m, M]) with m < a- < a+ < < P+ < M, and 

f{s)ds 



I{a_,a+,/3_,/3+) 



a+ ^/is- q;_)(s - q;+)(^- - s)(^+ - s) 
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Let a* and /3* be two reals such that m < a* < (3^: < M. Let e = (61,62) G with 
61 = — and 62 = /?+ — /?-. Then there exists a constant ( eM. such that 
/(«,)(! + 0(62)) log 61 + + 0(61)) log 62 



/3* - a* 

as a± — >■ a^: and (3± — )■ /3=k, i'o ?/za? 6 = (61, 62) — > (0+, 0" 



C + 0(61 log 61, 62 log 62), 



Proof. It follows by applying the first and third estimates of the previous lemma to the 
integrals JJ^ and J^' for some point 7 G («+, although before we must fix the lower 
limit of the first integral with the change x — a* = s — a+, and the upper limit of the second 
integral with the change x — /3* = s — /3_ . □ 

Lemma 25. Let KeCu^ = be a family of square linear systems defined for 6 > 0. 

(i) If the limits K = lime„^o+ and r = lime_^o+ exist, and K is nonsingular, then 

= w + 0(|K, -K|,|r, -r|), e ^ 0+, 

where u = K^^r is the unique solution of the nonsingular limit system Kw = r. 

(ii) If in addition, the matrix and the vector are differentiable at e = 0, then the 
solution uje also is differentiable at e = 0. To be more precise, if 

K, = K + eL + o(e), r, = r + eC + o(e), 6^0+, 
for some square matrix L and some vector (, then 

Ue = CO + eK + o(e), 6 — )• 0^, 
where u = K.^^t and k, = K.^^{( — Lw). 

Proof Both results follow directly from classical error bounds in numerical linear algebra. 
See, for instance, BM §2.7]. □ 

Appendix A.2. Properties of the rotation number 

Let us write the rotation number as the quotient p(A) = A(A)/2_ft'(A), where 
A(A) = / -=, K{X) 







andT(s) = {X - s){b - s){a - s). 



The study of the limit A — 0+ is easy. From lemma 22 we get the estimate A (A) 



2{aby^/^\^/^ + 0(A^/2), as A ^ 0+, whereas from lemma 20 we get that 



K{X)= +0(A), A^0+. 

Jb ^y s[s — b){a — s) 

By combining both estimates we get that p(A) = k^X^^'^ + 0(A'^/^), so limA_s.o+ p(A) = 0. 
Next, we consider A — 5^. After some computations based on lemma [23| we get 

K{b-e)= -c"^/2log6 + r/ + 0(6log6), 6^0+, 
A(6-e)= -c"^/^log6 + /i + 0(6log6), 6^0+, 
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where c = a — b, r] = fi + r], fj. = jl + fj., with f) = c log 4c, /i = c log 46, and 
1 1 \ ds 2 /"^ dx log4 



^Ja- s yja-h J s-h y/c Jq x + a/c a/c 
~ _ /"Y ^ 1 \ ds _ _ 2 /-^ dx _ 1 4c 

^ io - s ^Ja-h) b - s VcJ^ x + ^c ^c (i/a + yc)^ ' 

We have used the change x^ = a — s in both integrals. Let r^^, = c^/^r/ = log 16c and 
= c^/^/i = log 166c(a^/^ + c^/^)^^. Then we have the estimate 

K[b — e) 1 — c^/^r/log e + 0(e) 1 — r;* log e 



as e — 



0+. Thus, = (t], - /i,)/2 = log {{a/by/' + (c/6)i/2) = log {d + {<f - ifl^) 



acosh c?, where c? = {a/by/"^. This implies that cosh^ = a/b. Besides, estimate (A.l ) is the 
key to prove that the caustic parameter A° is exponentially close to b. Once fixed p° < 1/2, 
let A° G -E be the unique caustic parameter such that p(A° ) = p°, < e = 6 — A° ^ 1, and 



6 = log e. By finding 6 ^ = log e in estimate (A.l I, we get 



loge = 1/5 = r/. + ^^^ + 0(e). 

Using that /f^ = (r^* — p*)/2 and r^* = log 16c, we check that A° = 6 — e, with 

e = e^/^ = e''*-2«''/(i-2p°)+o{e) ^ lece"''^/^^/^"''") + O J^g-^''''/^^/^"'"")) , 

as p° — 7- (1/2)^. The limit A — t- 6+ is completely analogous. We omit the computations. 

With respect to the limit A — )■ we note that T{s) = {b — s){a — s)^ + 0(a — A) 
uniformly in the interval [0, b]. Hence, 

ds 2 I b 

A(A) = / ; h 0(a - A) = ^ atan \ - + 0(a-A), A a . 

Jo [a-sjy/b-s y/a-b \l a - b 

we get the estimate K{\) = n(a — 6)^^/^ + 0(a — A), as A — a^. 



21 



Besides, from lemma 

Therefore, p(A) = g + 0(a — A), as A — a~, where the limit value g G (0, 1/2) is defined by 
tan^ ng = b/{a — b). That is, sin^ ng = b/a. 

Appendix A. 3. Another characterization of the frequency 

We associate a "frequency" tu = zu{c) G M" to any vector c = (ci, . . . , C2„+i) G M^"+^ such 
that Co := < ci < ■ ■ ■ < C2„+i in the following way. First, we consider: 

• The polynomial T(s) = Yl^2i^{ci — s) G ]R2^^+l[■s], which is positive in the n + l intervals 
of the form (c2j, C2j+i); 

• The n + 1 linear functionals P{s) ^ ICj[P{s)] = {T{s))-^/^P{s)ds; 

• The n + 1 column vectors Kj = (/Cj[l], /Cj[s], . . . , /Cj[s"-i])* G M"; 

• The n X n nonsingular matrix K = {—Ki, . . . , (— 1)"_K'„); and 

• The linear functionals /C[P(s); to] = /Co[P(s)] + 2 Ej=i(-l)-'^^j''Cj [P(s)], for w G M". 



The frequency map for billiards inside ellipsoids 



42 



The hypothesis ci > is not essential to get a nonsingular matrix K, but it suffices to assume 
the strict inequalities Ci < ■ • ■ < C2„+i; see [l25l §111.3]. 



Lemma 26. There exists an unique G M" such that 

]C[P{s)- uj] = 0, VP(s) G Rn-i[s], (A.2) 

or equivalently, such that Kq + 2Ku; = 0, which is the matricial form of the linear system 
given in ([^. 



Proof. By taking the basis {1, s, . . . , s'^"^} of ]R„_i[s], we see that condition (A.2 1 is 
equivalent to the linear system Kq + 2X0; = 0. □ 
Therefore, condition ( A.2| ) is an equivalent characterization of the frequency. From now 



on, u = vj{c) stands for the frequency computed through the previous steps. 
Appendix A.4. Geodesic flow limit: c\ — t- 0^ 

Let Kg = 2(nr=2 c.)-'/'. Kg = {Kg, 0, . . . , 0) G M", and T^{s) = -s U^^^'ic, - s). 
Let K*^ be the n x n nonsingular matrix associated to the vector c*^ = (0, C2, . . . , C2„+i). Let 
g j^n ^Yyq unique solution of the linear system Kg + 2K.^k'^ = 0. Then 

Lo = K^cy^ + 0{cl^'^), ci^0+. (A.3) 

The proof is short. First, we note that T = T*^ + 0(ci) uniformly in [0, C2„+i]. Thus, 



using lemma 20 we get that K = K"^ + 0(ci) as ci — )■ 0+. And using lemma 22 we see that 
Kq = Kgcl'^ 0(ci/^) as ci 0+. Therefore, the linear systems Kg + 2K^k'^ = and 
Ci^^'^Kq + 2K(c7^^^tj) = are 0(ci)-close, being K*^ nonsingular, so (A.3) follows from 
the first item in lemma [25i 

Appendix A.5. Simple regular collapse: C21+1, C21 c* for some I = 1, . . . ,n 

Set = (ci,...,C2i_i,C2i+2,---,C2„+i) G M2n-i_ LetT^(s) = ni^2z,2«+i(ci - s) be 
the polynomial associated to c^. Let JCf and /C^ be the functionals associated to c^. Let 
00^ = (cjf , . . . , cj^) G M", where := (wf , . . . , ufi^.ujf^^, • • • , w^) = ^{c^) e M""^ is 
the frequency associated to c^, and ojf' G M is determined by 

1 + 2> / ^ / + , = 0- (A.4) 



^/W{s) j^,J.,, \c*-s\^/W{s) ^f=f^] 
Let e = C2/+1 - C21. Then 

w = + O (e) , C2«+i , C21 -^c*. (A.5) 



In order to prove this claim, we first observe that characterization (A.2) is equivalent to 
the system of n linear equations 

/C[(c* - s)s';u] = Ofori = 0,. . . ,n - 2 
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because {1, c* — s, . . . , (c* — ^} is abasis of ]R„_i[s]. Now, using lemmas 20 and 21 
deduce the estimates 

/Cf[sHO(e) ifj</, 
0(e) ifj=/, 
-/Cf_i[/] + 0(e) ifj>/; 
and 

-1/2 

/C,[l] 



we 



r.[(c* - s)s'] = r''+' fe!^|l£d^ 



7r(-T^(c*))""' + 0(e 



C2j 



\c* - s\^TR{s) 



) if J = I, 

+ 0(e) otherwise. 



Therefore, the hnear system (A. 6 1 is 0(e) -close to the nonsingular linear system 



IC^[s'; u^i] = for i = 0, . . . , n - 2 



condition (A. 4 1 



and the asymptotic formula (A. 5 1 follows from the first item in lemma 25 



Appendix A.6. Simple singular collapse: C21-1, C21 — )■ c* for some I = 1, . . . ,n 



Set c 



s 



[Ci, . . . , C2/_2, C21+1, 



. Let T^is 



n 



CO 



s) be 



the polynomial associated to c . Let JCj and /C be the functional associated to c . Let 

- (wf , . . . , ) G W\ where u^, := fwf -, , w,^, , wf) = w(c^) G M'^"^ and 



^. Let /C}^ and /C^ 
r\ where wj^ := (wf , . . . , a;f_i, wf^^, . . . , ) = w{c^) E 

1/2 if/ = l 
wfl]^ otherwise 

Let e = C2/ — C2i-i > and 5 = \ loge|^^ > 0. Then there exists E M" such that 



CO, 



CO = LO^ + 5k^ + o{S), 



C2l-l,C2l C . 



(A.7) 



To prove this claim, we set d = ^JT^{c*) > 0. We know that characterization (A. 2 1 is 
equivalent to the system of n linear equations 

IC[6d;io] = 
IC[{c* -s)s';uj] 



for z = 0, 



,n-2 



(A.8) 



because {Sd, c* — s, . . . , (c* — s)s" } is a basis of ]Rn_i[s]. Now, using lemmas 20 and 23 



we deduce the following asymptotic estimates. On the one hand, there exist some constants 
Co, Ci, e M such that 

'l + 05 + O(e) ifj = /-l,/. 



ICj[6d] = 6d]Cj[l] 



Cj6 + 0(e5) 



otherwise. 
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/C,[(c* 



S]S 



C2j+\ 



[c* — s)s^As 



/C.^kl + 0(e) 



Ms 



C2j \c*-s\^TS {s)+0{e) 



C2I-2 ^TS{s) 
s'ds 



In particular, /Ci_i[(c* - s)s*] - ICi[{c* - 



_ rc2i+i 



Jc2l-2 y/TS(s) 



+ 0(e) 

+ 0(6) 



0( 



if J < / - 1, 
if J =1-1, 

if J > /. 



/Cf_iM + 0(e) 



We assume now that / 7^ 1. The case / = 1 is studied later on. Since e <^ S, the linear 
system (A. 8 1 is 0(5)-close to the nonsingular linear system 

2(-l)'-H^f-i 







which in its turn is equivalent to the linear system 



C21+1 



Ms 



for i = 0, 



i-i 



K,^[s^-uj%] = Ofori = 0, 



,n-2 



(A.9) 



i-i- 



whose unique solution is uj^i = w{c°) and uj^ 

Thus, the asymptotic formula u = + 0{5) follows from the first item in lemma 25 



In fact, this result can be improved using the second item in lemma 25 It suffices to note that 
the linear system (A. 8 1 is not only 0(5) -equivalent to (A.9), but it is differentiable at 5 = 0. 
Hence, ( A.7 1 holds for some vector that could be explicitly computed in terms of the limit 
system and the constants Co, • • • , Cn- 

If / = 1, the linear system (A. 8 1 is 0(5) -equivalent to the nonsingular linear system 



1/2 



for 2 = n - 2 



and the proof ends with just the same arguments that for 1^1. We omit the details. 



Appendix A.7. Total regular collapse: C2Z+1, C21 — )■ for all I = 1, . . . ,n 

Let us study the case of n simultaneous collapses, all of them regular. That is, once fixed 
a vector c* = (c*,...,c*) G M" such that < ci < < ■ • • < c* , we study the 
asymptotic behavior of the frequency tu = zu(c) when C21+1, C21 — )• for all / = 1, . . . , n. 
Let CO = {coi, . . . , oJn) G M" be the vector whose components verify that < < 1/2 and 
sin^ TTUi = ci/c*i. Let e = (ei, . . . , e„) G with e; = C2/+1 — C21. Then 

io = u + 0{e), e-> (0+,...,0+). (A.IO) 

Let Qi = - ci Yli^i \c* - c;| > 0. Let {Pi{s), . . .,Pn{s)} be the basis of M„_i[s] 
univocally determined by the interpolating conditions 




Qi if j = I, 
otherwise. 
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That is, Pi{s) = (— 1)'"^a/c^ — ci Hi^/lc* — s). Using lemma 21 we get the estimates 



^^[m] = r + O(^) 1 = 2(-l)'-^ atan . /-^ + 0(e 







/C;[P/(s)] = TT + 0(e), and /Cj[Pi(s)] = 0(e) for j 7^ 0, /. Thus, the n x n linear system 

/C[Pj(s);w] = Ofor/ = l,...,n 
is 0(e) -close to the nonsingular decoupled linear system 

2(-l)'"^ l^atan ^ - ixuj^ = for / = 1, . . . , n, 

whose unique solution is given by tan^TrcD^ = ci/(c^ — ci), and so, by sin^ vro;; = ci/c^. 
Hence, the asymptotic formula ( A.lOp follows from the first item in lemma 25 



Appendix A.8. Total singular collapse: €21^1, C21 — )■ c* for all I = 1, . . . ,n 

Let Q = {1/2,..., 1/2) e W\ e = (ei, . . . , e„) G Wl, and 6 = {61, ...,6n) e W^, where 
Q = C21 - C2i_i and 61 = \ loge^l^^ Then 

CO = 10 + 0(6), e (0+,...,0+). (A.ll) 

Remark 14. By applying repeatedly the result on simple singular collapses, we see that 

lim ( ■ ■ ■ lim ( lim u] ] = to. 

In fact, these repeated limits can be taken in any order. Nevertheless, this result is weaker than 
estimate ( A.l 1| ), so we need a formal proof of the estimate. 



We consider the constants Qi = A/c2n+i — Hi^/ ~ c^l- Let {Pi(s), . . . , Puis)} be 
the basis of Mn-i [s] univocally determined by 



Qi if j = I, 
otherwise. 

Now, using corollary [24} we get that there exists some constants Qi E M such that 

l + Cii6i + o{6) ifj = l, 
Cji^i + o(^) otherwise, 
where < j < n and 1 < I < n. Therefore, the n x n linear system 

Si}C[Pi{s);uj] = Ofor / = l,...,n 
is 0(5) -close to the nonsingular linear system 
1 - 2cDi = 

2(-iy-\L0i-i - uji) = Ofor / = 2,...,n 



whose unique solution is cD = (1/2, ... , 1/2). Thus, the asymptotic formula (A.ll I follows 
from the first item in lemma |25| 
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Remark 15. The vectorial estimate (A.ll I can be refined in several ways. For instance, one 
can get the componentwise estimates wi = 1/2 + 0{5i) and ui = cui^i + 0(5/) for / > 1. In 
particular, coi = l/2 + 0(Si, . . . , Si). Even more, there exists anxn constant lower triangular 
matrix L such that 

co = LO + L5 + o{5), e ^ (0+,...,0+). 

We omit the proof, since we do not need this result and the computations are cumbersome. 

Appendix A.9. Asymptotic behavior of the function 

The function : (0, c) U (c, 6) ^ M verifies that /(A) + J{\)p^{\) + K{\)v^{\) = 0, where 
the coefficients /, J, K : (0, c) U (c, 6) — )■ M were given by 

= r ? — VTwrv ^(^^ = f -( — VTwrv ^^^^ = twt' 

Jq {a-s)y/T^{s) Jrn [a- s)^/Ty^[s) y/-T^{a) 

with Tx(s) = (A — s)(c — s){b — s), m = min(A,c), and m = max(A,c). Here, 
Px(A) = p(A; c, h) is the rotation function of the ellipse obtained by sectioning the ellipsoid 
Q with the coordinate plane {x = 0}. The asymptotic properties of rotation functions of 



billiards inside ellipses were established in proposition 10 

First, let us consider the case e := A — )■ 0+. Using lemmas 20 and 22 we get: 

• /(e) = + 0(e3/2), ^^^^^ ^ 2a-\hc)-^'^; 

• J(e) = Jo + 0(e), where Jq = -2 /^^(a - s)-\s{s - c){b - s))-^/'^<ls; 

• K{e) =Ko + 0(e), where Kq = 27r(a(a - c)(a - b))-^/^; 

• Px(e) = /t'^e^/2 + 0(e'^/2), where k,'^ = K^{b, c) can be found in proposition 10, and 

• zy^(e) = -(/o + JoK^)Koh^/^ + 0(e3/2) = 0{e^/^). It is possible to check that 
(Jo + Jok'^)Kq^ < 0, but we do not need it. 

Next, let us consider the case e := 6 — A — )■ 0^. We begin by computing the integral 

/•^ 2 

r(/3, a) := / . ry, = . ^ atan V/3/(a -/?), 

Jo {a-s)y/p-s \/a-p 

for any < /3 < «. Then it is immediate to check that 

ds r(/3,a_) -r(/3,a+) 



/o (a+ — s)(q;_ — s)v^/3 — s — a_ ' 

for any < /3 < a_ < a+. We also need the formula r(/3,a) = 27r(a — /3)~^/2f)(/3, a), 
where f)(/3, a) := lim^_^Q,- p(7; /3, a) is one of the limits of the rotation number described in 



proposition [10} Using these formulae, jointly with lemmas [20| and [21] we see that: 

• /(6 _ e) = + 0(e), where U = 27r(a - b)-^ {{b - c)-^/^g{c, b) - {a - c)-^/^g{c, a)) ; 

• J{b -e) = J, + 0(e), where J, = -27r(a - b)-\b - c)-^/^; 

• K{b -e) = + 0(e), where = 2n{a - b)-\a - c)-^l'^\ 

• p^{b - e) = p,{b) + 0(e) = p(6; c, b) + 0(e) = f?(c, b) + 0(e); and 
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Vy,{h - e) = ^(c, a) + 0(e) = p(a; c, a) + 0(e) = py{a) + 0(e). 

The estimates in the limit e := m — m — )• 0+, which means A — )■ c, are: 

/(c ± e) = —(a — c)"^(6 — c)~^/^ loge + /i + O(eloge), where /x is a constant that can 
be exactly computed from lemma 23 , 



• J(c ± e) = — (a — c) ^(6 — c) log e + r/ + 0(e log e), where r/ is a constant that can 
be exactly computed from lemma 23 , 

• K{c ± e) = 27r(a - c)-i(a - h)-^'^ + 0(e); 

• px(c ± e) = 1/2 + K-^log'-^e + 0(log~^e), where = K^{c,b) = acosh(6/c)^/^ 
according to proposition [lOl and 

• z/,(c±e) = (a -6)1/2 ((^ _ c)(ry _ ;^) - 2(6 - c)-!/^^^) /27r + 0(log"^ e). After some 
tedious, but simple, computations, one gets that i^x(c) = g{b, a) = p{a; 6, a) = Pz{a). 

Appendix A. 1 0. Asymptotic behavior of the function z/y 

The function z/y : (6, a) R verifies that /(A) + J(A)py(A) + K(A)z/y(A) = 0, where the 
coefficients I, J, K : {b,a) R were given by 



(6-3)^1^' 7a (3-6)^1^^' ^/=%(by 

with Ty(s) = (c — s)(A — s)(a — s). Here, Py(A) = p(A; c, a) is the rotation function of the 
ellipse obtained by sectioning the ellipsoid Q with the coordinate plane {y = 0}. 
We begin with the limit e := A — 6 — )• 0+. Using lemmas [20| and 23 we see that: 



• /(6 + e) = Io + 0(e), where Jq = J^ib - s)-^/\c - s)-'/\a - sy^^^ds; 

• J{b + e) = 27r(a - b)-'/\b - c)-VVi/2 + 0(1); 

• K{b + e) = -27r(a - b)-^/\b - c)-!/^!/^; 

• py(6 + e) = py(6) + 0(e); and 

• z/y(6 + e) = p^(c) + 0(e^/2)^ ^j^ce py(6) = p(6; c, a) = p(c; 6, a) = p^(c). 

Next, let us consider the case e := a — A — )■ 0^, which is similar to the limit 
lime_>o+ '^x(6 — e) studied in the previous subsection, so we need the same simple integrals. 
Using them, jointly with lemmas 20 and|2T] we get: 



/(a - e) = + 0(e), where h = 27r(a - 6)-^ ((6 - c)-^^^q{c, 6) - (a - c)-^/'^g{c, a)) ; 
J{a - e) = + 0(e), where = 2Tx{a - 6)"^(a - c)-^/^; 
K{a -e) = K, + 0(e), where = -27v{a - b)-\b - c)-^/^; 
Py(a — e) = g{c, a) + 0(e); and 

z/y(a - e) = g{c, 6) + 0(e) = p(6; c, 6) + 0(e) = p,{b) + 0(e). 
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AppendixA.il. Asymptotic behavior of the function 
The function : (0, c) U (c, 6) — )■ M verifies that 

/(A) + J(A)p,(m) + K{X)u,{X) = 0, (A.12) 
where the coefficients J^K : (0, c) U (c, 6) — )■ M were given by 

with Tz(s) = (m — — s){a — s), m = min(A,c), and fn = max(A,c). Here, 
Pz(A) = p(A; b, a) is the rotation function of the ellipse obtained by sectioning the ellipsoid Q 
with the coordinate plane {z = 0}. 

First, let us consider the case e := A — )■ 0+. Using lemmas 20 and 22 we see that: 

/(e) = /oei/2 + 0(e3/2), ^^^^^ ^ 2c'\ah)-^/^- 

J(e) = Jo + 0(e), where Jq = 2 /^"(s - cY^{s{s - b){a - s))-^/^ds; 

K{e) =Ko + 0(e), where Kq = -27r(c(6 - c)(a - c))-!/^; 

Pz(m) = pz(min(e, c)) = pz(e) = K'^t^l'^ + 0(e^/2), where the constant k*^ = k'^(&, a) 
can be found in proposition 10, and 

z/,(e) = -(/o + Jo«:^)^o"'eV2 + = 0(eV2), ^j^^ (Jq + Jo«:^)^o"' < 0- 

The estimates in the limit e := m — m — )• 0+, which means A — c, are: 
/(c ± e) = 7r(a - c)-i/2(^ _ ^ o(l); see lemma 

J(c±e) = 0(1); 

ir(c ± e) = -27r(a - c)-'/^{b - c)-'/\-'/^ + 0{e'/^y, 

Pzim) = Pz(min(c ± e, c)) = Pz(c) + 0(e), since Pz(A) is analytic at A = c; and 
z/,(c±e) = l/2 + 0(ei/2). 

Next, we consider the case e := 6 — A — > 0+. Using lemmas 20 and 21 we get: 
/(6-e) = 0(l); 

J{b - e) = J,e-i/2 + 0(1), where J, = 2n{a - b)-^/\b - c)-^/^; 
K{b - e) = K*e-i/2 + 0(eV2), where iT, = -27r(a - 6)-V2(5 _ c)-'/^; 
Pzim) = Pz(min(6 - e, c)) = Pz(c); and 
z/z(6-e)=Pz(c) + 0(eV2). 
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Appendix B. A topological lemma 

We recall that the complement of any Jordan curve X in the plane has two distinct 
connected components. One of them is bounded and simply connected (the interior, denoted 
by Bx) and the other is unbounded (the exterior, denoted by Ux)- 
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Lemma 27. Let X and Y be two Jordan curves o/M^. If f : Bx — ?■ is a bounded local 
homeomorphism that has a continuous extension to the boundary X such that f{X) C Y, 
then f : Bx — > By is a global homeomorphism. 

Proof. We note that W = f{Bx) is a non-empty open bounded subset of such that 

dW = df{Bx) C fidBx) = fix) c Y. 

Next, we are going to prove that W = By- Using that dW C Y, we deduce that the 
intersection WD By (respectively, WCiUy) is open and closed in By (respectively, in Uy), so 
it is either the empty set or the whole interior (respectively, exterior). Therefore, we deduce 
that: 1) n = 0, because W is bounded; 2) ly fl F = 0, because W is open; and 3) 
W n By = By, because W is open and non-empty. That is, f{Bx) = W = By. 

Once we know that / : Bx — ^ By is a surjective local homeomorphism, we deduce 
from covering space theory that it is a global homeomorphism. It suffices to realize that Bx 
is connected and open, and By is simply connected. □ 

In particular, if / : Bx — ?• is smooth or analytic, then its inverse is also smooth or 
analytic. This means that if / is a local diffeomorphism whose image is bounded and that 
has a continuous extension to the boundary X such that /(X) C Y, then / : Bx By is a 
global diffeomorphism. 
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